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ABSTRACT 


A flight control system for use in air-to-air combat simula- 
tion has been designed. The input to the flight control system are 
commanded bank angle and angle of attack, the output are commands 
to the control surface actuators such that the commanded values 
will be achieved in near minimum time and sideslip is controlled to 
remai n smal 1 . 

For the longitudinal direction, a conventional linear control 
system with gains scheduled as a function of dynamic pressure is 
employed. For the lateral direction, a novel control system, con- 
sisting of a linear portion for small bank angle errors and a bang- 
bang control system for large errors and error rates is employed. 


PREFACE 


The work reported herein was performed under Contract NASl- 
13773 between January 1975 and January 1976. The study was con- 
ducted under the direction of W. Hankins, Simulation and Human 
Factors Branch, Langley Research Center, NASA. 

This study was performed at Decision Science, Inc. with 
George H. Burgin acting as responsible project manager and David 
M. Egglestonasprincipal investigator. 

The purpose of this contract was twofold. The primary goal 
was to design and to demonstrate a control system that would ac- 
cept, from the "Adaptive flaneuverinq Logic" (AfIL) program commands 
for angle of attack and for bank angle and translate these com- 
mands into control surface deflections. 

The secondary goal of this contract was to make available to 
the Simulation and Human Factors Branch at LRC a package of com- 
puter programs which can be used in the design of such a control 
system for other aircraft than F-4. 

Both these goals have been achieved. The aircraft driven by 
the AML program and controlled by the control system described in 
this report not only performs as well in simulated air combat than 
the previous performance model, but in some instances oerforms 
superior to the performance model. 

The computer programs described in this report demonstrate 
that the secondary goal has also been achieved. Enough documen- 
tation of these programs is provided in this report so that pros- 
pective users of these programs should have no problems in using 
them to design similar control systems for different aircraft. 
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DESIGN OF AN ALL-ATTITUDE FLIGHT CONTROL SYSTEM TO 
EXECUTE COMMANDED BANK ANGLES AND ANGLES OF ATTACK 

By George H. Burgin and David M. Eggleston 
Decision Science, Inc. 

SUMMARY 

This final report on Contract NASl-13773 is organized in three 
different major parts plus an Appendix. 

In the first part, the general concepts and the philosophy 
that were used in developing a flight control system for the F-4, 
as represented by the equations of motion programmed on the DMS, 
are described. The flight control system accepts as input com- 
manded angle of attack and bank angle and determines from these 
and the current state variables of the controlled aircraft the 
required stabilator, aileron and rudder deflections. The general 
six degree of freedom equations of motion are presented, as well 
as the linearized equations for the longitudinal and the lateral 
axes. Discussed are also difficulties of general nature associ- 
ated with the roll maneuver. 

The second part specializes in the above concepts for an 
F-4, with special emphasis on the adverse yaw and inertial 
coupling. An attempt is made to formulate a set of performance 
specifications for such a control system, and the associated dif- 
ficulties are discussed. Finally, the control system as designed 
for longitudinal and lateral axis is given and the method of de- 
termining its gains is explained and some typical response time 
histories are presented. 

The third part of the report is devoted to a description of 
the computer programs developed under this contract. Detailed 
flow charts of the major programs and subroutines are presented, 
the input requi reinents for each individual program are listed and 
instructions for the use of these programs are presented. Two 
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major programs and a set of five supporting programs were develop- 
ed. The program ATDYN is a versatile, general purpose program for 
linear control system transfer function analysis and the program 
AML75 employs the developed control system to control an F-4 
driven by the AML program in an air-to-air combat simulation. The 
supporting five programs are used for pole and zero determination 
in closed-loop root locus studies for determining reasonable 
values for control system gains. 
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INTRODUCTION 


It is a well known fact that piloting an aircraft in close-in 
air-to-air combat is one of the most demanding tasks for a pilot. 

A thorough understanding of the complex dynamic relationships 
governing the flight in air-to-air combat is required for pilots 
to be trained effectively for this difficult task. This same in- 
sight into flight dynamics is a prerequisite for the desigh of new 
fighter aircraft, associated weapons systems and flight control 
systems . 

Considerable progress has been made during the last years in 
providing the aircraft design engineers and the pilots with air- 
to-air combat simulations. These range from relatively unsophis- 
ticated flight simulators to air-combat simulation, whefe real 
aircraft, equipped with range instrumentation, "^fy over* ah air- 
combat maneuvering range. Flight instructors can monitor the en- 
gagement from ground consoles; the engagement may also be replayed 
for debriefing purposes. 

The Differential Maneuvering Simulator (DMS) at the NASA, 
Langley Research Center (LRC) is a good example of a flight simu- 
lator which can be used for air-to-air combat simulation and which 
provides a very high degree of realism. It is used in LRC's 
Advanced Fighter Technology program and is described in Reference 
1. Among the research tools available on that simulator is a 
computer driven, interactive opponent, selecting its maneuvers 
according to an Adaptive Maneuvering Logic (AML) (References 2, 3 
and 4). This AML driven opponent, often called the "iron pilot", 
is indeed a worthy opponent and usually "beats" its human adver- 
sary. 

However, the mathematical model used in the AML program is a 
performance model. That means, thrust, lift, drag and turning 
capabilities are represented at their maximum, steady state values. 
The AML driven aircraft therefore does not adequately represent the 
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transient behavior and the associated handling qualities of the 
aircraft. As a consequence, the AML driven aircraft never becomes 
uncontrollable or unstable. This gives the "iron pilot" a definite 
advantage, most pronounced in the low speed, low dynamic pressure 
flight regimes. 

The objective of this study was to design a flight control 
system which can be used to replace the performance model of the 
present AML program by equations of motion identical to those rep- 
resenting the simulated aircraft flown by the pilot. Thus, the 
AML driven aircraft will have identical handling qualities and per- 
formance as the piloted aircraft. It will no longer be a perfectly 
stable aircraft, but it will have the same spin and departure char- 
At tM tircrift fltw by tht iMMiafi pilot. 

TM thrtt pristry control vtriabltf In air coabat are bank 

angle, load factor and ffiagnittrde of thrust. Thrust Is simulated 
by assuming either idle, military or afterburner thrust. Bank 
angle and load factor have to be controlled by commanding input to 
the actuators of the control surfaces of the aircraft, typically to 

stabllotor, alloron, iptlltr ati4 rudder actuators. 

SMMiffr tilt mi pnoffM •n«cet*» a tactical decision, typi- 
cally once every second, it determines a desired maneuverpl ane ro- 
tation angle and a desired load factor. These commands can be 
translated into a desired bank angle arid desired angle of attack. 
?he^ prdbl em Laddresaed in ithi sssftudy can be sumfiiarized as follows: 

For commanded bank angles (({) ) and commanded 
angles of attack (a^,), which are given at one 
second intervals, find a control law that 
will make the aircraft follow these commands 
in some "optimal" manner. 

Keeping in mind that such an attitude control system has to 
operate over practically the entire flight envelope of a given 
high performance fighter aircraft, this seems like a formidable 
task. However, one must realize that the problem is not quite so 
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horrendous as it might seem at first glance for the following 
reasons : 

(a) The dynamics of the plant (here the aircraft) are accurately 
known. Unlike in the design of a control system for a real 
aircraft, where stability and control derivatives are only 
estimated and have uncertainties associated with them, in the 
case considered here all these coefficients are exactly known. 
Also known is the exact form of the differential equations 
describing the motion of the aircraft to be controlled. 

(b) The aircraft is represented as a rigid body. 

(c) All aircraft state variables are accessible without noise, 
without measurement errors and without sensor dynamics. 

(d ) The number of nonlinearities in the system is relatively 
small. 

This report is organized in such a way that it proceeds from 
the general topics, applicable for the control system design for 
any fighter aircraft to results which are specific for the F-4 and 
finally to the detailed descriptions of the computer programs de- 
veloped in this study. 

To design an efficient flight control system, the basic, open 
loop aircraft dynamics must first be undeirstdo'd thoroughly. The 
open loop characteristics of the aircraft suggest certain require- 
ments of the control system. For instance, strong allverse yaw may 
suggest certain strategies for ai 1 eron- rudder control for rolling 
under high angles of attack. 

The first part of the report explains the basic philosophy 
and the procedure used in the design of the flight control system. 
The results presented in this part are aircraft independent and 
they can be used for the design of a similar system for some other 
advanced fighter aircraft. 

The second part presents results specifically for the F-4. 
First, the linearized open and closed-loop behavior of the F-4 
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aircraft is explained, followed by results for the full set of the 
nonlinear equations of motion. 

Also included in this second part is a discussion of a prelim- 
inary set of performance specifications. In the course of the work 
performed under this contract, it became quite clear that formu- 
lating a set of performance specifications for an aircraft/control 
system for such a complex mission as air-to-air combat becomes a 
formidable task which in itself would require more resources than 
were allocated to the entire contract. 

The present tendency in aircraft specifications is still to 
describe the aircraft as an equivalent second order system. It is 
left to the airframe designer and to the control system designer 
to develop an aircraft which resembles a second order systems. 

To the authors of this report, such an approach appears unaccept- 
able for air-to-air combat maneuvering. Not even the most gifted 
and ingenious designer will ever build an aircraft, which, when 
performing the drastic maneuvers required in air-to-air combat, 
will behave like a second order system! 

The third part of the report is devoted to the description of 
the computer programs developed under this contract. Flow charts 
and detailed input and output descriptions are provided. In ad- 
dition to the two major programs, AML75 and ATDYN, a set of five 
small , indiyiiMal programs for linear control system design are 
described. 

Of the two major programs, AML75 may be considered as the end 
product of this contract. It allows exercise of the developed 
control system in any one of the three following modes: 

(1) Single aircraft, driven by an inputted command 
sequence of angle of attack and bank angle at 
one second intervals. 

(2) AML dri ven target aircraft against canned 
trajectory attacker. 
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(3) AML driven target aircraft against AML driven 
attacker using the old performance math model. 

The program ATDYN is a tool for linear transfer function an- 
alysis. It incorporates the subroutine TRANS, a versatile and ve 
powerful transfer function analysis program. This program was de 
veloped and used for many years at General Dynamics/Astronautics 
and later General Dynami cs/Convai r for control system design and 
analysis of the Atlas and Centaur missiles. It determines poles 
and zeros of closed-loop systems without requiring that the cor- 
responding equations be set up in first order form. 
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SYMBOLS, 


SIGN CONVENTIONS AND TERMINOLOGY 


Since the equations of motion as used in this study are iden- 
tical to the equations of motion as documented in Reference 5, the 
symbols used in this> report are, wherever possible, compatible with 
the ones used in Reference 5. For the sake bf'completeness, most 
of these symbol s are ' 1 i sted agai n in the following list: 


^long 


Coefficient matrix for linearized longitudinal 
equations . 


1at 


Coefficient matrix for linearized lateral 
equati ons . 


b 

^long 
^lat, I 
®lat, 2 



Wing span. 

Vector of forcing function coefficients for 
longitudinal equations (fig). 

Vector of forcing function coefficients for 
lateral equations (6^). 

a 

Vector of forcing function coefficients for 
lateral equation (fip). 

Mean aerodynamic cord. 

Drag coefficient. 

Lift coefficient. 

V. . 

Rolling-, pitching*, and yawing-moment co- 
efficients , . resp^ecti vely . 

Longitudinal-, lateral-<i and verti cal -force 
coefficients, respectively. 

Total external force vector. 

X, y and z components in body axis system of 
the vector F. 
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g 




^yy’ ^zz 


xz 


la 


‘Pa 


K 


‘pe 

p<i> 


‘Ra 

‘Re 


■^R(|) 

thru 


L 

L, M, N 


M 

m 


Acceleration due to gravity. 

External moment vector. 

Angular momentum vector. 

Components of vector h. 

A1 ti tude. 

Moments of inertia about body x, y and z axes, 
respecti vely . 

Product of inertia about body x and z axes. 
Gain for integral a-feedback. 

Gain for integral (j>-feedback. 

Gain for proportional a-feedback. 

Gain for proportional 3-feedback. 

Gain for proportional (f>-feedback . 

Gain for rate a-feedback. 

Gain for rate 3-feedback. 

Gain for rate <|)-feedback. 

Coefficients in linearized equations of mo- 
tion (1 thru 19 longitudinal, 30 to 51 for 
1 ateral ) . 

Lift. 

Rolling, pitching and yawing moment, respec- 
ti vely. 

Distance between center of gravity and aero- 
dynamic center of tail. 

Mach number. 

Total aircraft mass. 
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P, q. r 


q 

s 

s 

T 

u, V, w 


V 

X. Y, Z 



a 

“( 

а, 

3 

б . 


w 


’SB 

*SP 


Rotational rates about x-body, y-body and 
z-body axes, respectively. 

Dynamic pressure. 

Wing area. 

Complex Laplace transform variable. 

Thrust . 

Components of velocity along x-body, y-body 
and z-body axes, respectively. (In sections 
about linear analysis, u, v and w denote small 
changes of these variables from a steady- 
state reference value.) 

Velocity vector of aircraft center of gravity. 

Aerodynamic forces along the x-body, y-body 
and z-body axes, respectively. 

Aerodynamic forces along the x, y and z axes 
in the stability axes coordinate system. 

Aircraft x and y coordinates in earth fixed 
coordinate system. 

Body angle of attack. 

Commanded body angle of attack. 

Wing angle of attack. 

Angle of attack error. 

Sideslip angle. 

Aileron deflection. 

Rudder deflection. 

Stabilator deflection 

Speed brake deflection. 

Spoiler deflect! on . 
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t|>, e, (|) 

<I>E 

on 

■^lat 


Ai r dens i ty . 

Euler yaw, pitch and roll angle, respectively. 
Commanded Euler roll angle 
Error between actual and commanded angle <(». 
Stabilator actuator lag. 

Aileron actuator lag. 


w Mrenft wtiflit. 

At In fttftrtftct S, tilt tquttloiif In tk<t r«^ert yst d1«tn- 

si.onal control derivatives and static force and moment derivatives. 
For example. 


C 


L 


3 


63 


Nondi mensi onal dynamic derivatives are used, such as 

C„ _ 2V 

F“ 6p • 


Sign convention for control surface deflections. 


positive for trailing edge of left aileron down, 


6 positive for right spoiler up. 

5 p 


Ai 1 eron : 

Spoi 1 er : 

Rudder: 6^ positive for trailing edge left. 

Stabilator: 6^ positive for trailing edge down. 


11 


A note on terminology .- Throughtout this report, the expres- 
sion "old equations of motion" refers to the equations of motion as 
they are used in the present (1975) version of the AML program on 
the DMS. Characteristic for these equations of motion is the fact, 
that no moments of the aircraft are calculated. The attitude of 
the aircraft is obtained in the following way: Forces acting on 

the center of gravity are calculated under the assumption that the 
aircraft has a given angle of attack and zero sideslip angle. By 
integration of these forces, the velocity vector at the end of an 
integration step is obtained. The requirement of the net forces 
lying in a given plane, a given angle of attack and a vanishing 
sideslip angle then define uniquely the attitude that the aircraft 
should have at the end of this integration cycle. Taking the dif- 
ference between the actual and the desired aircraft attitude de- 
fines the desired attitude change. Postulating that this attitude 
change can be acconplished by three rotati ons about the aircraft body 
axes during one integration stepsize defines the desired body ro- 
tational rates p, q and r. These rotational rates are then pro- 
cessed through a digital filter which provides the actual body ro- 
tational rates . These rates are integrated to provide the actual 
attitude at the end of an integration step. 

The expression "new equations of motion" refers to the com- 
bination of the nonlinear, six degrees of freedom equations of 
motion as used in the aircraft simulation on the DMS and the dif- 
ferential equations of the flight control system as developed in 
this study. The new equations of motion can be considered as 
describing the (nonlinear^ transfer function between the input 
commands (angle of attack, bank angle and thrust level) and the 
aircraft motion. 
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PART I: SENERAL CONCEPTS 
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CONTROL SYSTEM DESIGN PHILOSOPHY AND PROCESS 


The task of executing AML commanded maneuvers is clearly a 
problem of an automatic control system synthesis. Several syn- 
thesis methods have evolved over the last years, one of the best 
summaries of these methods can be found in a paper by Whitaker 
and Potter (Ref. 6). Although this paper is almost ten years old, 
it is the opinion of the authors of this report that most of the 
facts and opinions expressed in this paper are still true today. 

The syiitNsis methods may be classified into two broad cate- 
gories. The older procedure, which may be considered to be an 
outgrowth of classical feedback theory, i s an iterative trial and 
error process. It operates mainly in the frequency domain, ad- 
justs gains, time constants, etc. until the location of the roots 
of the transfer function of the closed-loop system indicate proper 
frequencies and dampings. It is interesting to note that the MIL- 
SPECS MIL-F-8785, "Military Specifications - Flying Qualities of 
Piloted Airplanes" rely to a great extent on these same concepts 
of natural frequency and damping. This trial and error procedure 
may, today, seem somewhat unscientific, but when executed by an 
experienced design engineer, it generally takes only a few i ter- 
ations to obtain desirable root locations. 

The newer approach to control system design is based on what 
is generally called modern control theory. Here, the design of a 
control system is started by defining a performance index and the 
control system is then designed such that the performance index 
is minimized. For linear, second order systems, analytical methods 
are known to minimize certain performance indices. But even for 
this simple design task, the computational effort may become 
uneconomical ly fapge. In ; addi ti on , an ai rcraf t , especial- 
ly when flown in air-to-air combat, cannot be represented as a 
linear second order system or as several, independent second order 
systems. The nonlinear effects of the cross coupling terms in the 
equations become very significant. In a later section of this 


report, the problems arising from inertial cross coupling are ex- 
plained in more detail. 

The fact that the actual equations of motion of the 
aircraft are nonlinear, of course, would not exclude the use of 
the principle of minimizing a certain performance index; it would 
only exclude the possibility of finding an analytical solution to 
the problem. The almost unsurmountabl e difficulty appears to lie 
in the definition of meaningful performance indices. All the clas- 
ical methods of linearization are based upon perturbations about a 
nominal this nominal flight con- 

di tion is 'trimm^^ Teyel , unaccelerated fl i|ht. Such a flight 
condition may prevail just prior to the beginning of an engagement, 
but thereafter, during the dogfight, it is unlikely that the 
straight and level flight is ever resumed. Typical air combat 
maneuvers may find the aircraft initially in a steep dive, banked, 
say, 90 degrees to the left, and a relatively low angle of attack. 
The required maneuver may then be a 90-degree right roll combined 
with a pitch to an angle of attack which corresponds to a maximum 
lift coefficient. It is conceivable that for this particular ma- 
neuver, one might find a meaningful and realistic performance in- 
dex and then iterate on the design of a control system for that 
particular flight situation. The next maneuver, which may be 
commanded, may be a change in bank angle without a change in angle 
of attack. It is unlikely that the same performance criterion, as 
in the previous maneuver, would be desirable in this situation. 

This illustrates the difficulties associated with specifying per- 
formance criteria. Since these should be made a function of flight 
condition, may be dynamic pressure, due to the nonlinearity of the 
equatlj^s^ of. they would also have to be a function of the 

command.input, single axis commands yetsus multiple axes command 
and small changes 'versus large changes. 

Adaptive Control Systems .- Adaptive Control Systems are char- 
acterized by incorporating a mechanism which automatically adjusts 
the control system parameters (usually gains) to adapt the control 
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system to either changes in the external environment or the flight 
conditon or to changes in the aircraft dynamics. Ever since the 
first adaptive control systems were used with the X-15 aircraft, 
they have been quite popular with experimental aircraft. The main 
advantage of an Adaptive Control System is that the aircraft dy- 
namics (that is, its stability and control derivatives) do not have 
to be known before the control system is actually used in the air- 
craft. In the particular task here, however, the dynamics of the 
aircraft is exactly known, since it is given by the DMS equations 
of motion. Gonsequently. there is no need to have an Adaptive 
Control Systent. ^ 

Design Process .- The fastest way to obtain a working control 
system capable of executing command sequences generated by the AML 
program appeared to be to use classical, linear control systems 
resembling classical pitch and yaw dampers in stability augmenta- 
tion systems. Figure 1 shows, in the form of a flow chart, the 
design process. The design process starts with selecting repre- 
sentative flight conditions. Three conditions were selected, all 
at 15,000 feet straight and level flight; a very low dynamic pres- 
sure (q=150 psi), a medium dynamic pressure (q=300 psi), and a 
high dynamic pressure (q=840 psi). One might also consider ad- 
ditional flight conditions with low, medium and high dynamic pres- 
sure, one at close to sea level and another one at real high al- 
titude, say, 35,000 feet. To familiarize one's self with the basic 
properties of the open loop airframe at these flight conditions, 
it is recommended that one obtain the transfer functions a(s)/6s(s) 
and 4)(s)/6aCs) as well as the location of the roots of these trans- 
fer functions. This provides a first indication of what one might 
reasonably expect in terms of pitch and roll response. 

Next, a basic longitudinal and lateral control system is 
selected and initial gain values are determined by engineering 
judgment, by comparison with similar control systems and by in- 
spection of the open loop response. Ideally, then, root loci plots 
are developed by varying one gain at a time. If it is possible to 
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Figure 1 (Concluded) 




find gains which will result in a satisfactory distribution of the 
roots, the design proceeds to the next step. If, however, the 
roots cannot be moved to desirable locations, the control system 
itself, rather than just its gains, need modification. Either more 
variables should be fed back or additional elements, like lead/lag 
networks should be added and new gains determined. 

After a feedback system has been found whose roots appear to 
satisfy the desired response characteristic, a crucial step follows, 
and that is to determine how well the nonlinear equations of motion 
can be controlled with the control system obtained by the linear 
analysis. For small perturbations, the nonlinear equations of 
motion should behave very closely to the linear case unless cer- 
tain types of nonlinearities, such as deadbands and hysteresis are 
present. If such nonlinearities exist, it is recommended to first 
remove those from the nonlinear program to ol»taiN a naaningful 
coiiparitgfi llntar and nonllntar analysis. 

If, for SMall aispTItutfn, singit axis disturbances, the re- 
sponse of the nonlinear system is drastically different from the 
reponse of the linear system, the most probable cause is an error 
in one or both of the two programs. Careful comparison of the 
computer program against the original equations and comparison of 
all the coefficients between the two programs may reveal the error. 
Other potential problem areas exist in the mathematical formula- 
tion of the differential equations of the control system and in an 
inadequate method of integration or an unsuitable integration 
stepsize. On the other hand, since the programs for the linear 
analysis contain several critical portions, the following checks 
should be made: Are the matrices, from which the eigenvalues are 

determined, set up properly? Are the signs of all the feedback 
paths correct? Do the eigenvalues found by the computer program 
really satisfy the equation [A] - X[I] = 0? Note that the com- 
parison between linear and nonlinear equations offers a valuable 
tool for program checkout and debugging. 
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Once the response of the linear and nonlinear systems agree 
well enough so that it seems likely that the two programs are error 
free, one may proceed to the next step which basically investigates 
how much are linear and nonlinear systems different for: 

- small amplitude longitudinal and lateral 
command input combined. 

- large amplitude single axis command input. 

- large amplitude longitudinal and lateral 
command input combined. 

Each one of these three tests may require slight modifications of 
the gains or, in extreme cases, even the control system. 

After a satisfactory response for single command input (lat- 
eral and longitudinal simultaneously and individually) has been 
obtained, the program is exercised by command input sequences. The 
objective here is to investigate how the system behaves if it ob- 
tains a new command input before the response to the previous command 
input has settled to steady state. During this process, the air- 
craft is also operated at flight conditions different from the ones 
assumed at the beginning of the design process. Thus, unsatisfac- 
tory response characteristics may be revealed at certain flight 
conditions. The closed loop response at these critical flight 
conditions may then be investigated by repeating the entire pro- 
cess outlined so far, starting with determining all the coeffic- 
ients of the aircraft by the program ATDYN. 

If the control system has passed all the tests so far, it is 
reasonable to assume that it has now the capability to "fly" the 
aircraft with a command sequence actually generated by the AML. A 
canned trajectory is read in for the attacker aircraft, the AML 
program controls the target. If a reasonable attacker trajectory 
is selected, the target aircraft has to fly under flight conditions 
typical for air combat maneuvering. By this time, the main ob- 
jective as far as the control system is concerned, should be to 
prevent the aircraft from departing. Note that this is now the 
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combined res pons ibillty of the AML program and the control system. 
The AML program should command maneuvers such that the aircraft 
will not enter flight conditions with extremely low dynamic pres- 
sure in which the aircraft cannot be controlled. The control sys- 
tem should at all times execute the AML commanded maneuvers, if 
physically possible, in such a way as to keep the aircraft under 
control. Specifically, simultaneous roll and pitch maneuvers must 
be executed so that good use is made of the specific aerodynamic 
properties of a particular aircraft. 

The final step in the design process is to use the AML pro- 
gram with the new equations of motion against an AML program with 
the old equations of motion. The primary objective of this mode 
is to compare performance between the new and old eouations of 
motion, while secondary objectives are to obtain estimates about 
computer time required for implementation of the new eouations of 
motion in real time en ^the DMS and to obtain additional information 
concerning the performance of the AML program with the- new equa-' " 
tions of motion. Most important here is again the question of 
stability. Any time the AML-driven aircraft should become unstable, 
an investigation about its cause must be made. The crucial ques- 
tion obviously is: Should the AML program have avoided getting 
into an undesirable flight condition or was the control system not 
capabile of performing a legitimate AML maneuver in a stable manner? 

Figure 2 illustrates the individual programs v/hich support the 
linear analysis and gain selection. 

During the work performed under this contract, the situation 
that the AML-driven aircraft, running either against canned tra- 
jectories or against another AML-driven aircraft, became unstable, 
did never occur. However, only a very limited number of such runs 
were made, and it is possible that in a more extensive checkout, 
occasional departures of the AML-driven aircraft might occur. 
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Figure 2. - USE OF LINEAR ANALYSIS AND CONTROL PROGRAMS 
STEPS IN ANALYSIS AND GAIN SELECTION FOR 
A SPECIFIC FLIGHT CONDITION 
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EQUATIONS OF MOTION 


The general, six degrees of freedom, equations 
of motion for a rigid body 

The most general form of the rigid body equations, expressed 
in a space-fixed coordinate system, are the well known two vector 
equati ons 


F = ^ (mv) 


G = 


dh 

^[t 


where 


F 


V 

h 

m 


resultant external force vector acting on center of 
mass 

resultant external moment vector acting about center 
of mass 

velocity vector of center of mass 

angular momentum vector 

mass 


Transforming these equations from a space-fixed reference frame 
to one rotati ng with the aircraft, the classical Euler equations of 
motion are obtained, in vector form: 

F = m^+mwxv 


fih 

G = + w X h 


or, expressed in the six corresponding scalar equations; 
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= m(0 + qw 

- rv ) 


= m('5 + ru 

- pw) 

fz 

= m(i^ + pv 

- q'u) 

L 

= * '<!'z 


M 


- ph. 

N 

• 

= h + ph 

z ^ y 



Expressing the angular momentum components and their derivatives in 
terms of moments of inertia* we obtain the standard body-axis 
rigid-body equations of motion as used on the DMS (and on most 
other flight simulators): 

F 

U = — + rv - qw 

m ^ 

V = pw - ru 

w * ^ + qu - pv 

p = L ^ ^ ^ , qr <'yy.~ 

XX XX XX 

q = « + (^2 . p2) ^ . pr i' . ^LX l . 'zz) 

yy yy yy 

N . /I ^xz __(^yy ” ^xx) 

r = ^ — + (p - qr) j pq ■■ t - 

zz zz zz 
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The above form of the equations of motion is most useful for 
analysis; for programming on a digital computer, they have to be 
rearranged because the equation for p contains f and the equation 
for r contains p. This computational difficulty can be overcome by 
substituting the r equation into the p equation and then solving 
for p , which yields : 

P =[l - pq +(n . pq(ly^-l„) - qr} 



The equation for r needs no reformulation because, when r is cal- 
culated, p is already known. 

Linearized Equations for Longitudinal Direction 

The linearized equations of aircraft motion used in this study 
are derived in the book by Blakelock, Ref. 7 . The equations des- 
cribing longitudinal (pitch) motions given on page 21 of Ref. 7 
may be written in the form (for definition of axes, see Figure 3) 

U + Kg U + a + a Kg9 + Kg 6 = (x-axis force) 

• ' ' . 

Kg U + Kg a + K^q a + K^^0 + K-| 2 ® “ *^13^s (z-axis force) 

*^14 ^ *^15 “ *^16 *^17^ ^18® ^ *^19^s (pitch moment) 


with the coefficients K^ through K^g being defined as: 


K 


1 


mVy 

sq 


K 3 = - 2 ^ (neglected) 

'^T 





sq au 
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Kr ~ “T\ 7 ~ ^v_ (nsglected) 


K,, = V- C = C 


^ COS0 = -C., cos 9 

® Sq w 


" "^m 


K? C (neglected) 


‘15 


•2Vt 


C (u+Au) -C (u) 

''s ' -<:z = ^ 

u Au 


‘16 " ‘"'m 


a 


~ ^1 0 — ^7 » ^7 ^ ~ — 7 “ 

^ * *T^JL ^ Jt I' c 

T & & IT & Sqc 


“^lO ■ 

a a 


18 2V^ 


*^11 " ■'^1 " 2 Vt ^z 

• q 


“^IQ “ ^m 


K ^2 = ^ sine 
Sq 


where 


• u * 

u = 
U. = 


£ 

U 

x-axis velocity change 

x-axis velocity at steady flight condition 


“ = w <■“> 


'a = angle of attack perturbation 
0 = pitch angle perturbation 
= stabilator deflection 

m = total aircraft mass 
Vj = total velocity 

S = wing planform reference area 
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— 12 
q = dynamic pressure = j pVj 

X = aircraft forward body axis 
y = aircraft right-wing body axis 
z = aircraft downward body axis 
"c = mean aerodynamic chord 
V = Y-axis velocity change 
w = x-axis velocity change 

tj = length between c.g. and aerodynamic center of tail 

The coefficients in the K array are defined in terms of the 
nondi mensi onal stability derivatives such as, for example, 

1 ^^x 

= •?— - — . Since an adequate discussion of stability deriva- 

Sq “ 

tives for aircraft is beyond the scope of this report it is as- 
sumed that readers lacking familiarity with them can refer to 
Ref. 7. 

The term K 2 involves the partial derivative which depends 

upon both drag and thrust variations with forward speed. Since air 

combat usually requires either full afterburner or idle thrust, the 

corresponding variation in thrust with forward speed was evaluated 

and included in x for either afterburner or idle thrust in AML75. 

3u 3F 

In CMDSEQ the value of — ^ is calculated based upon afterburner 
thrust. 


Linearized Equations for Lateral Direction 

The linearized equations for the 1 ateral -di recti onal axes 
(roll-yaw) are given on page 116 of Blakelock's book as 

*^30*^ *^31*^^ * ’^ 32 '^ ’^ 33 '^ *^34^ “ *^35'^a *^36^R 

K37(j> + K^gij) + Kggil^ + *^43^R 

'^44‘f* '^46'^ *^47’*' *^48^ *^49^ ' '^50'^a '^51^R 

force) 
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where 


^30 


Sqb 
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'32 


^33 


2Vt 

'xz 


Sqb 

b 


•7V7 % 




*^35 " ^-1 


*^36 " 


^37 *^32 


K 
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2V. 


'39 


2z 

Sqb 


'40 




Ixx” x-axis moment ..of tnerti a 

b = wtng span 
6. = aileron defjleGtfbrv 

O 

6|^ = rxidder deflecti 0 :n . , 


*^41 ■ '^n, 


•^42 ■ ^n 


^43 = 


R 


•^44 = rv^ Sp 


K 


45 = C = coS(|) cos(|) 
y<f, Sq 

= -K^ coS(j) 


'46 


•^1 - ^ s , 


'^47 “ ' ■’'12 


'^48 ■ *^1 


•^49 ' S, 


'<50 ' Sj, 


Kci ' C 

51 y{ 


^ 2 z~ moment of inertia 

g = standard gravitational ac- 
cel erati on 

I = y_z product of inertia = 

Jr ^ 
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THE ROLL MANEUVER 


Maneuver Limiting Factors 

To roll the aircraft fast and accurately is one of the key 
tasks in air combat maneuvering. The roll maneuver is convention- 
ally defined as a rotation of the aircraft about its longitudinal 
axis with an angular velocity p. As W. H. Phillips has shown (ref. 
11), even smal 1 disturbances about a condition of rapid, steady 
rolling may cause longitudinal and lateral instabilities. If 
there is a significant angle or angular rate between the velocity 
vector and the aircraft longitudinal axis, the situation gets even 
more complicated. The following two situations are quite impor- 
tant for the design of a control system and are, therefore, dis- 
cussed in some detail in the following: 1) Initiation of a roll 

from a flight condition with a large angle of attack, and 2) Sim- 
ultaneous rotation about the x-body and the y-body axis (simul- 
taneous roll and pitch). 


R«11 NtMiivtr Under Nifh Anflt of Attack 

If an aircraft performed a pure rotation about its longi- 
tudinal axis only, after a rotation of ninety degrees the initial 
angle of attack would be pure sideslip and after 180 degrees, the 
initially positive angle of attack would be a negative angle of 
attack with the same magnitude as the initial angle of attack. To 
verify the above statement , consider an a ire raft flying initially 
with wings level, pointing north, and an angle of attack of a . 

Its direction cosine matrix D then is 


cosa 

0 

" s i n a 

0 

1 

0 

si not 

0 

cosa 


The velocity vector in inertial coordinates is: 
32 


Si nee 


it follows that 


/ K 


I 


and by definition 


a 


0 

0 


[D] 


atan (^) 
a s i n ( ■ 


Xg cosa \ 
0 


Xg sina 


'Total 


-) 


a 

3 


a 

0 


Let now the aircraft perform a rotation of 90 degrees about 
its x-body axis, then 


[D] 


and consequently 

a = 
8 


cosa 
si na 
0 


0 

a 


0 

0 

-1 


-sina 

cosa 

0 


The importance of this phenomenon, as far as control system 
design is concerned, lies in the fact that rolling under angle of 
attack causes sideslip even if no side forces were generated, 
simply by the geometrical relationships. Rolling under angles of 
attack therefore requires simultaneous yawing. To keep a roll 
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maneuver, which is initiated with no angle of attack, coordinated, 
a yawing moment must be generated to compensate for the yaw moment 
due to C^p (yaw due to roll) and due to ^ and (adverse 

yaw due to aileron and yaw due to spoilers). Small rudder deflec- 
tions will be sufficient to keep the net yaw moment close to zero. 

For large angles of attack, the situation is quite different. 
Not only is the adverse yaw stronger (requiring more rudder de- 
flection), but as it was shown above, additional yawing moment 
must be generated to yaw the aircraft in such a fashion that the 
angle of attack remains in the aircraft symmetry plane and does not 
degenerate into sideslip. 

In Section II of this report these relationships are illus- 
trated by some numerical examples for the F-4. 

Simultaneous Roll and Pitch Maneuver, Cross Coupling 

The above discussion assumed that during the roll maneuver, no 
change of the angle of attack was desired, that is, that the angu- 
lar rate about the body y-axis remained essentially zero. Often, 
in air-to-air combat, the optimum maneuver consists of changing 
both the bank angle and the angle of attack. During such a maneu- 
ver, both roll rate and pitch rate are nonzero and, in order to 
keep the maneuver coordinated, the yaw rate will also be nonzero. 
Inspection of the equations of motion shows that under this con- 
dition, the terms qr ( ^yy ~ ^zz ) (etc) become effective. This 

^ XX 

effect is generally called inertial cross-coupling. Special at- 
tention has been paid in the past to the q equation. The inertial 
cross coupling term there is 


<'xx - 'zz) 

‘yy 

Since in most modern fighters, I is considerably smaller than 

A A 

I^ 2 » high roll rates couple strongly into the longitudinal motion 
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and some aircraft become actually unstable at high roll rates. 

But even the cross coupling term in the p equation, that is 

qr (I - I ) 
yy 22^ 


I 


XX 


can become significant in the design of a control system for air 
combat maneuvering. Even when and are of the same order, 
since they both are large, the difference between the two may be 
quite significant, and since it is divided by the small moment of 
inertia about the longitudinal axis, its effect is amplified. Once 
the existence of a significant cross-coupling term in the foil ac- 
celeration equation for a given aircraft has been established, the 
control system designer has to take some corrective action. 

For instance, for an F- 4 , where is greater than I , the 
cross coupling term for the roll acceleration equation is of the 
form 

-K qr 

where -K represents the negative quantity 


The sign of the cross-coupling term obviously changes with the sign 
of q. If an increase in angle of attack is desired, q will be pos- 
itive, for a decrease of the angle of attack, q will be negative. 

It seems reasonable to assume that, for a desired roll moment, 
depending on the sign of q, the cross-coupling term miaht add to or 
oppose the desired roll moment. 


One might therefore consider to incorporate some logic into 
the lateral and longitudinal control systems which would cause the 
aircraft to roll and pitch simultaneously if the cross-coupling 
term assists the roll moment but perform roll and pitch sequential- 
ly if the cross-coupling term is of opposite sign to the desired 
roll moment. 
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SPECIFIC CHARACTERISTICS 
FOR THE F-4 
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OPEN LOOP CHARACTERISTICS OF THE F-4 


General Characteristics 

To effectively synthesize an attitude control system, it is 
necessary to know and understand the dynamic and aerodynamic pro- 
perties of the basic aircraft. Of prime importance are those af- 
fecting fast pitch and roll maneuvering over the entire flight en- 
velope encountered in air-to-air combat. 

The original plan was to present time histories for control 
surface step input, first stabilator, aileron and rudder steps 
alone and then combined step input. While the open loop poles of 
the linearized transfer functions indicate what type of response 
may be expected, it would be interesting to compare the response of 
the full set of nonlinear equations of motion with the results ex- 
pected from the poles. Several such comparisons were made, but, as 
pointed out In the preface, the limited resources did not permit ad- 
equate documentation of these results. While they were important and 
useful to arrive at the final control system, once the control sys- 
tem is designed, it is no longer necessary to have access to these 
intermediate results. 

As just an example. Figures 4 and 5 show typical open loop 
responses for an aircraft with an initial sideslip angle and all 
control surfaces in their neutral trim position, both responses at 
15,000 feet altitude, the first with M=0.32 and the second with 
M=0.8. Additional open loop characteristics can be obtained from 
Reference 5. 


Adverse Yaw 

At the beginning of this study, it was believed that due to 
the extremely high adverse yaw of the F-4 roll commands at high 
angles of attack could not be initiated by aileron control. How- 
ever, inspection of the pertinent control derivatives at Mach num- 
bers below 1.1 shows that the following is true; 
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Figure 5. - Open Loop Response at M=0.8 at 15,000 Feet 
Sideslip Angle = 5° at t=0. 
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(a) Cp is positive tof all flight regimes. 

(b) Cp is almost independent of the angle of attack. 

^6a 

(c) Cp is positive for all flight regimes. 

^^6sp 

(d) Cp decreases strongly with increasing angle of attack. 

^^6sp 

As a consequence, even under high angles of attack, aileron 
(and the associated spoiler) deflection will create a rolling moment 
in the desired direction. Note that the rudder deflection occurs 
automatically due to the B-feedback loop into the rudder. Fortun- 
ately, the control loops as designed are capable of properly com- 
pensating the adverse yaw effect. 
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SOME PERFORMANCE SPECIFICATIONS 


Problems Associated with Performance Specifications 

Specifying performance characteristics for an ai rcraf t/control 
system configuration may be compared to specifying handling quali- 
ties of an aircraft. The well known specification, MlL-F-8785 
"Military Specification - Flying Qualities of Piloted Airplanes" 
(see, for instance. Ref, 8), illustrates a possible approach to a 
set of specifications. 

The main difficulty in defining specifications for aircraft 
perform^ince during air combat maneuvering is the fact that the 
aircraft does not behave like a second order linear system and 
that, therefore, response characteristics are dependent on the am- 
plitude of the maneuver under consideration. Not only are, of 
course, the basic equations of motion nonlinear, but severe non- 
linearities exist in the limits of control surface deflections and 
possibly in the rates of these deflections. 

While, for instance, in a linear system, the response for a 
commanded one-degree change of angle of attack would be the ident- 
ical time function as for a two-degree change in angle of attack, 
with simply a different scale on the ordinate axis, this is not 
the case for a nonlinear system. 

A first approach would, therefore, be to set up performance 
specifications for realtively small commanded changes. As long as 
certain types of nonlinearities, such as dead zones and hysterisis, 
are absent, these "small command" specifications certainly give 
some indication about the aircraft performance. 

In air-to-air combat, however, the commanded changes are not 
small, but they may be extremely large, like from trim angle of 
attack to maximum angle of attack. It might, therefore, be more 
meaningful to specify certain constraints which should be satis- 
fied in such extreme maneuvers. Even if such specifications were 
set up, say maximum rise time, etc., to demonstrate that these 
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specifications are satisfied would require substantial effort (max. 
angle of attack for each flight condition has to be established first). 

For the lateral direction, "small commands" would be of the 
order of a few degrees, while again the maneuvers commanded in air- 
to-air combat often much larger. Time to roll from level flight to 
say 90 degrees might be a reasonable specification here. 

Note that in all the above discussion, single axis command 
inputs were assumed. Clearly, the performance criteria set up in 
this single axis mode have to be relaxed if simultaneous command 
input are considered. The amount of individual specifications re- 
quired becomes nov/ quite large, because already in the single axis 
mode, the performance criteria are at least a function of dynamic 
pressure and of angle of attack. 

It is suggested that the entire subject of performance spec- 
ification be studied in a separate effort, not necessarily assoc- 
iated with the AML program. Special emphasis should be given to 
performance specifications at high angles of attack, where roll 
reversal, reduced roll control power and, particularly on F-4 air- 
craft, spin susceptibility may limit the maneuvering capabilities. 

For an attitude control system responding to commands from the 
AML program, the most important requirement is fast response to 
achieve the commanded maneuvers. The system must also be stable 
for continuous flow of input in all combinations of input magni- 
tude and frequency. 

It is quite difficult to define performance specifications for 
an aircraft control system over the entire flight envelope of air 
combat. The aircraft dynamics are complex and the responses are 
nonlinear. Even at the same speed and altitude, the dynamics and 
control characteristics still depend upon angle of attack, pitch 
angle and roll angle. The time required for a roll maneuver is a 
nonlinear function of the magnitude and sign of the desired roll 
maneuver. 
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At very low speeds, which can occur as the result of air com- 
bat maneuvering, roll control of the F-4 becomes very difficult. 

The adverse yaw of the ailerons at low speed and high angles of 
attack is so powerful that even hard rudder is insufficient to 
maintain control during attitude maneuvers. 

To completely characterize and specify the controlled respon- 
ses of the aircraft over this entire flight envelope in a restric- 
tive but achievable manner is a task worthy of the entire resources 
of this study. 

In order to reduce the specifications to manageable form the 
following restrictions were imposed: 

(1) Speed range - from Mach 0.3 to Mach 1.1. 

(2) Specifications are given only for the aircraft in 
trimmed, level flight. 

(3) Specifications are given only for a limited set 
of flight conditions. 

Comments on restriction 1 .- While the control systems devel- 
oped will successfully fly an F-4 slower than a velocity corres- 
_ 2 

ponding to q = 150 Ib/ft a control system optimized for fast, 
effective attitude control at very low q would require a different 
approach to aileron and rudder control. This would be an inter- 
esting study but is beyond the scope of the present resources. 

Comments on restrictions 2 and 3 .- The aircraft dynamics vary 
significantly with angle of attack and even change somewhat with 
pitch and roll angle. The computer time available under this con- 
tract was not sufficient to allow simulation of the controlled 
responses for a large number of flight conditions after the best 
performing control system and gains had been found. This simula- 
tion would be necessary to evaluate the performance achieved for 
other flight conditions. 

The flight conditions chosen are representative of the entire 
flight envelope. Rise time and overshoot for step input commands 
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for other flight conditions nearby is expected to be similar. 

Longitudinal Rise Time for Step Input 

The most important specification for an air combat attitude 
control system is rise time for a commanded step input. Figure 6 
shows the spec|ft'Cation and the achieved values of rise time as a 
function of dynamic pressure at an altitude of 15,000 feet. Also 
shown for comparison is the computed period of the short period 
mode of the basic aircraft. Sequences of pitch command step input 
were givep whicjR Included a range of sizes of steps from a few 
degrees up to the maximum change from the trim angle of attack 
possible without overstressing the aircraft. The points enclosed 
by a square represent the average rise time while the highest and 
lowest recorded rise times are shown using horizontal bars above 
and below the squares. 

The observed rise times were usually less than the time of one 
period of the short period mode. The fastest rise times are a- 
chieved near Mach=.8. The control system should never exceed the 
specification line which is formed of straight line segments be- 
tween the points; 

q = 150 Ib/ft^ tr = 3.0 sec 

q = 850 Ib/ft^ tr = 1.6 sec 

q = 1500 Ib/ft^ tr =1.4 sec 

Figure 7 shows the level flight trim angle of attack and 
the maxii%um allowed angle of attack for flight of the F-4 at 15,000 
feet. This- i nformation is necessary to find the magnitude of al- 
lowed change in angle of attack. 

Lateral Rise Time for Step Input 

The rise time in roll for flight at 15,000 feet is shown as a 
function of dynamic pressure in Figure 8. Some points obtained 
in simulation runs with the bang-bang roll control system are also 
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Figure 6 . Allowable Rise Time for 
Angle of Attack Response 



Ib/ft 









shown. The bang-bang system gives extremely fast and effective 
rol 1 control . 

In Figure 9 is shown the observed rise time for the roll 
control system at Mach=1.0 at 15,000 feet as a function of the com- 
manded roll angle in degrees. The control system chooses the pro- 
per direction in which to roll so that the roll maneuver executed 
in one second never exceeds 180 degrees. The rise time appears 
nearly constant out to about 50 degrees which is due partly to the 
effect of the linear control system for small roll errors. With 
the bang-bang roll control system in effect the rise time increases 
nonlinearly as a function of commanded roll angle. 

Although specifications were not developed for flight at 
1000 feet and 35,000 feet the curves of wing trim angle of attack 
for the F-4 are shown in Figures 10 and 11 versus Mach number. 
These would be useful in establishing and verifying rise time 
specifications for pitch and roll step input at these flight con- 
dlti ons . 


Other Specifications for Controlled Response 

Although the rise time is critical in specifying control sys- 
tem response for air combat attitude maneuvering, other performance 
measures may be necessary or desirable. Primary among these are 
the overshoot, the settling time, and the steady-state error. 

The amount of overshoot is controlled mainly by the natural 
damping of the aircraft and the artificial damping due to pitch 
rate q, roll rate <{), and sideslip angle rate 3. Since we had es- 
sentially perfect, noise-free values of these rate signals, over- 
shoot in responses to step input could be controlled with no dif- 
ficulties whatever by raising the rate gain coefficients. Note, 
however, that raising rate gain in a loop results in increasing 
the rise time. 

With good rate signals, as in this study, the overshoot re- 
sulting from single-axis step input cbrnmands was usually less thati 
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Wing Angle of Attack 




Wing Angle of Attack 


tgu're 11 .» Wing Trim Angle of Attack 
for Level Flight 
F-4 at H=35,000 ft. 
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two degrees in pitch and five degrees in roiil , regardless of the 
magnitude of the input step. Thus, specification values given be- 
low should always be satisfied when the control gains are properly 
adjusted. 

Specifications for Single-Axis Step Input Response 

For step commands Aa^ in pitch: 

Overshoot must be less than five degrees in magnitude 
for step input command Aa less than 15 degrees to final 
values of a in the range of 0°<a^20° starting from trim- 
med level flight. 

For step commands ^ in roll: 

Overshoot must be less than 20 degrees in magnitude for 
step input A(J>^ up to 180 degrees starting from a condi- 
tion of trimmed level flight. 

For combined maneuvers simultaneously in pitch and roll the 
control system incorporates logic to avoid large roll errors due 
to rigid body coupling. Before this logic was implemented, roll 
overshoots as high as 360 degrees were observed during extreme 
roll rate maneuvers with combined pitch commanded changes. With 
the addition of the control logic and the development of the bang- 
bang control system for roll, the observed overshoots for combined 
maneuvers were dramatically reduced. Due to the limited resources 
of this study, no specifications were developed and verified for 
combined pitch-roll maneuvers. 

The settling time in response to a step input is usually de- 
fined as the time elapsed before the response enters and remains 
inside a region +5 percent of the input step magnitude centered 
on the steady state output. For the purposes of air-to-air combat 
it would be preferable to make the settling time within +^10 per- 
cent of the steady state value or to give the settling time in 
terms of absolute errors in pitch or roll. A pitch error of +1.5 
degrees and a roll error of +3 degrees could be recommended. 
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THE LONGITUDINAL CONTROL SYSTEM 


Longitudinal Control for Linearized Equations 

Since one of the ground rules for the control system design 
was the unlimited accessability to all state variables for feedback 
purposes, the most logical choice for selecting a feedback path was 
to feed back the very variable which was to be controlled, that is, 
the angle of attack. Note that the linearized equations for the 
longitudinal direction have been written in such a form that a is 
one of the state variables. 

The control system used for the longitudinal or pitch axis is 
shown in Figure 12. The proportional gain Kp^ represents the number 
of degrees the stabilator should be deflected per degree of angle 
of attack error a^. At low speeds and low dynamic pressure ^ a 
larger stabilator deflection is needed to exert fast control, while 
at high speed and large dynamic pressure, a considerably lower Kp^ 
must be used to avoid overstressing the aircraft. From a speed of 
response point of view Kp^ should be as high as the other contraints 
will allow. 

The integral gain is necessary to remove steady-state er- 
ror because Kp^ cannot be raised high enough to do so. If Kj^ is 
too low, steady state errors will persist indefinitely, while if it 
is too high, loop i ns tabi 1 i ty may result. 

Rate feedback gain Kp^ provides an error signal proportional 
to ^ which damps the pitch motion and provides lead compensation 
to the longitudinal control loop. If Kp^ is too high, pitch re- 
sponse will become very slow and sluggish. This gain should be 
set only as high as necessary to provide damping of pitch axis 
motions. 

The program. LONLOP (Longitudinal Loop) calculates the dynam- 
ics of the pitch control loop for any set of gains. 

The transfer function is negative since positive stabilator is 
defined as down elevator (stick forward) which results in a reduced 
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Figure 12. - Longitudinal Control Loop 





angle of attack. 

This transfer function shows a set of zeros that very nearly 
cancels the phugoid poles, indicating that the phugoid mode is not 
very important in angle of attack response to stabilator deflection. 
The phugoid mode is usually of very low frequency, in this case 

,/Wi th a period of ~ = 91 seconds. The period of the short 

period mode, on the other hand, is ^ = 3 seconds. If these 
roots are plotted in the complex plane, they appear as in Figure 13. 

A typical response of a system to a step input is shown below: 



U.* Typical Rttpoiitt to Stop ComoiiO Inpyt 


Here, t^ is the delay time, t^, is the rise time, and t^ is the set- 
,./.'tling:,t;1ra|i. ■ 

i For the attM^ delay plus* rise 

time must t»e as Is phySTd^l ly |>oss11>lfr,and the^ 

must stilt be kdpt Vies on ably smattvTKeU time 


57 


importance. Small variations in a about the proper mean value will 
have very little effect on the aircraft trajectory in air combat. 


Procedure for Determining Longitudinal Gains 

The first step in determining desirable gains is to run ATDYN 
for the flight condition of interest and obtain a data deck of the 
K array and the printout of roots. This output includes both lon- 
gitudinal and lateral data. 

The transfer function for ^ ^ ^ between stabilator and angle 

of attack is given in the printout section entitled LONG. WITH 
ALPHA/DELS ZEROS as the ratio of the third order polynomial to the 
fourth order polynomial. For the standard conditions (V=634.7 ft/ 
sec for Mach .6, H=15,000 ft, a=6=3.3375°) there results 

g(s) ^ (-.0773 - .04426s - 14.18s^ - .0957s^) 

^s^ (.0306 + .0224s + 6.446s^ + 2.0039s^ + 1 .2634s^) 


which in factored form using corresponding roots becomes 


a 

6 



-.0957 [(s+. 00154)^ + .0719^] (s+148.2) 



phugoid poles short-period poles 


Steady-state errors are important and cannot be allowed. They 
would significantly compromise the ability of the aircraft to fol- 
low AML commanded maneuvers. 

Once the transfer function is understood the next step is 
to solve for thi roots of the combined (control system + aircraft), 
or closed-loop, as. portrayed in Figure 13. This computation is done 
by LONLOP. Input data for LONLOP consists Of the K-array deck plus 
a set of longitudinal loop gains Kj^, The first run of 

LONLOP can be done with these gains set to zero. The resulting 
poles (denominator roots) should be the same as those obtained in 
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ATDYN, and this serves as a check. 

To obtain very fast response relatively high values of Kp^ and 
*^Ia necessary. Kp^ as high as four can be used at very low 
speed, while Kp^=2 may be quite high enough at high speed. Values 
of Kj^ usually run between 25% and 150% of the Kp^ used. The rate 
gain sbppld be as low as possible to result in_ keeping the 
damping of the closed loop short period poles at least as high as 
C=.25. 

It is quite helpful to plot the location of the closed-loop 
poles for ranges of one gain with the others held either constant 
or in a constafit ratio. In the present case, since the controller 

Ki Kd s + K. -Kd (s+fjL) 

transfer fiincti on is ^ = -( — ^ -) = “ 

s 

Kj 

it is convenient to keep the ratio -rr^ constant. 

This fixes the zero location due to this term at the point 

-Ki 

s = A plot of closed-loop root locations with this ratio 

Pa 

fixed and a series of values of Kp^then gives the root locus for 
proportional gain changes. Although not necessary for the use of 
these programs a knowledge of the procedures and meaning of the 
root locus design technique is definitely useful. Readers lacking 
this background can obtain a brief introduction from Appendix C of 
Reference 7 (Blakelock). 

Once a set of suitable gains is found they should be tested 
using the program version of NONLIN to find the response to a step 
change in a. Usually the damping obtained in these nonlinear model 
runs is higher than the linear program results indicate. 

Gains s.uitable at high speed will not be appropriate at low 
speed, so that for the complete flight envelope at least three and 
possibly more sets of gains, each applicable over some portion of 
the flight envelopb, will be necessary for longitudinal control of 
a gi ven a1 rcraft. 
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Longitudinal Control fon Nonlinear Equations 


Basically the same control system as described under the lin- 
ear longitudinal control is implemented in the AML75 program. 

Since angle of attack in the nonlinear equations is not one of the 
state variables, its rate is computed as 

uq-vp + g (t^+qS[C^ sinap+C^, cosap]/'^ + CD33}-uw/^^ 
u+q Sg C^. cosa /j^ 

The variable > = p+(q sin4» + r cos<(i) tane Is alSo computed 
for 1 ater use. 

A control logic was incorporated to avoid commanding extreme 
pitch maneuvers while relatively violent roll maneuvers were in 
progress. This logic is one possible way to avoid the pitch-roll 
coupling problem. 

The logic implemented first calculates the magnitude of the 
present roll errpr - <j) | and the magnitude of the roll rate er- 
ror |4» - <i> 1 Ul- If the conditions | c|) - 4> | >30 degrees and 

|(|)1 > 50 deg/sec are simultaneously satisfied LOCNTR goes into a 
"pitch slowly" mode such that 

“cMD = “ + 1 

This means that the pitch control aims for a point halfway 

between the existing angle of attack and the desired angle of 

attack a^. 

c 

An additional logic was ^added wliioh ensures that during such 
a relatiyely violent maneuver no .adglesi-o-f. at'^^^^^ greatet 
than Id degrees This logic rndmentar^kl^- brings 

aircraft to a low angle of attack region where roll is 

most effective. 
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When the error and error integral terms have been computed the 
stabilator deflection is calculated: 




Finally, the limits of stabilator motion are imposed, namely 
-23 deg < 6^ < + 7 deg. 

Scheduling of the gains as a function of dynamic pressure was 

implemented to provide fast pitch control over widely different 

flight conditions. For dynamic pressures 150 <^<520 Ib/ft^ the 
longitudinal gains read in as VAR(3), VAR(6), and VAR(IO) in the 
input data deck. are used. For convenience the longitudinal gains 
are written in the vector form (Kp^, Kj^, Kp^).w 

The gains for ^ < 150 Ib/ft^ are set at (3, 3, .6) and for 

q > 520 Ib/ft^ are set at (1, 2, .25) in LONCTR and can be changed 

only by altering the corresponding statements in the program. 
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THE LATERAL CONTROL SYSTEMS 


Linear Lateral Control 

The aircraft dynamics and consequently the control systems 
necessary for the lateral axes are more complex than for the lon- 
gitudinal axis. Since in this study all state variables are known, 
control of roll was based on knowing roll angle 4> and roll rate 
Control of rudder was based on knowing sideslip angle 3 and its 
rate B. 

rather complex nonlinear model of an aircraft the 
first approach was to determine what performance could be achieved 
with simple linear control and more particularly to find those con- 
ditions under which these basic linear controllers became marginal 

p|J ■ 1 

or exhibited poor performance or instability. 

The linear control systems developed for rudder and aileron 
control are shown in Figure 15. The 3 loop was set up originally 
without rate feedback. Large error excursions during combined ma- 
neuvers resulted, and the addition of 3 rate feedback brought al- 
most direct control of Dutch roll damping and greatly improved 
performance . 

Roll control was obtained using proportional, integral, and 
rate feedback. Integral control was required to reduce roll steady 

state error to zero. 

Quit# tarly In tills study it was found that a linear controller 

with proper gains would give very good roll control for snail step 
inputs, yet for large inputs with the same gains gave slow response, 
large overshoots, and generally poor control. Negative rate feed- 
back is also inappropriate for large roll maneuvers because it 
results in reducing the maximum roll rate and thus slowing down 
roll maneuvering. For these reasons, and because of the need to 
maneuver as fast as possible in air combat, a bang-bang control 
system was developed for aileron control. 
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PROP. AND RATE 

INTEGR. LI KIT 



Figure 15.- Linear Lateral Control Loops 









The Bang-Bang Control System 
for Aileron Control 

It is easily proven that for time-optimal control of a system 
with limited control authority the controls should always be de- 
flected to their stops* Thus, achieving a new roll orientation 
would involve applying full aileron in the direction to reduce 
roll error, followed by switching to full aileron in the opposite 
direction at the correct moment to result in zero error simulta-r 
neously with zero error rate. 

The general problem of time-optimal attitude maneuvering of 
an ai rcraf t f s qui te complex, and will not be treated here. It is 
apparent that a time-optimal attitude control system would have to 
give the fastest possible response in single-axis maneuvers as 
well as in combined pitch and roll maneuvers. Thus, time-optimal 
roll control would be an important segment of the problem. 

As a first approximation in developing the bang-bang control 
system, the assumption was made that the aircraft in its roll be- 
havior may b6 approximated as a simple double-integral plant. Note 
that by definition, the error is: 

<()e = <t>c - 4> 

and, consequently , the ordinate in the phase plane is 



Figure 16 a shows the phase plane plot for such a plant, for 
any given initial condition for roll error and error rate, these 
curves show the response to full negative aileron (solid lines) 
and to full positive aileron (broken lines). Mote that the curve 
contains the loci of all points for which, if full positive 
aileron is applied^ the system will simultaneously reach zero er- 
ror rate. Similarly, y_ is the location of all points with that 
same property for full negative aileron. 
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Assume that the system initially does not find itself on 
either the y+ the y_ curve, but anywhere within the region 
labeled R_ in Figure 16b. The basic idea of bang-bang control then, 
of course, is to apply full negative aileron until the solid line 
through our initial point has brought the system onto the curve 
the switching curve. At this moment, full positive aileron is 
applied until both error and error rate vanish simultaneously at 
which time aileron is neutralized. If the system finds itself at 
time zero somewhere in the region labelled R^, the same process 
with reversed signs is applied. 

Still under the assumption of a pure double integral plant, we 
can easily calculate the switching curves. This is done by back- 
tracking the switching curve from the origin; note that in the 
following, time is reversed (see Reference 9): 

.. _ '"max 


where L _ is the assumed rolling moment, constant during the con- 

nia X , 

trol maneuver. Since the "initial" condition i s <|)g = 0 and = 0, 
it follows that: 

ie(t) - t 

^ XX 


and 


♦ e(t) 


‘max 


I 


XX 



Eliminating t, we obtain 




e 



^e 
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Note that, if the three following assumptions were true, we 
would have a perfect time-optimal control system: 

1. Roll dynamics is pure double integration plant. 

2. Roll moment remains constant during the maneuver under 
i nvesti gati on . 

3. Aileron can be "switched" from full positive to full 
negative and reverse instantaneously. 

Since the three assumptions listed above are, of course, not 
satisfied in reality, we have to investigate how this affects the 
control system. Jhe main effect of assumptions one and two not 
being satisfied is a change in the shape of the phase plane curves. 

The response of a fighter aircraft to hard right aileron ex- 
hibits high roll acceleration as long as the roll rate is small, 
but soon approaches a maximum roll rate as a result of natural 
roll damping from the wing and limited aileron effectiveness. 

This effect is shown in the phase plane plot of Figure 17. 

If the aileron is instantaneously reversed when the aircraft 
is rolling at maximum roll rate, the roll deceleration is very 
large because both the aileron and the natural roll damping are 
then acting in the same direction. This is shown in the trajec- 
tory ABC, where at point B the aileron is switched from 6, = + 30 

a 

degrees to 6^ = -30 degrees. 

a 

Theoretically, to achieve time optimal roll control, the 
aileron would have to be reversed at exactly the right moment. 

The point of switching in the t))^, phase plane depends on the 
roll acceleration that will be attained after switching. Under 
the ground rules of this study (unlimited access to all state 
variables and complete knowledge of all the stability and control 
derivatives), one can compute the exact roll acceleration (decel- 
eration) obtained from full opposite aileron at any instant. The 
switching logic implemented here assumes this roll acceleration to 
remain constant during the time of switching to the time where the 
phase plane trajectory would pass through the origin. Section B-C 
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Area in which Linear 
Control Operates 



Figure 17.- Phase Plane for Aileron Control 




of trajectory in Figure 17. 

Note that in reality, the roll moment does not remain constant 
after the switching, but the roll deceleration will become smaller 
as the roll rate decreases. That means, the trajectory will not 
actually follow BC, but, say, BD. Fortunately, this is of little 
consequence, because at point D, the control system quits using 
bang-bang control and the linear control system takes over. 

One might ask the question, what happens if the trajectory, 
after the switching point B, misses the linear region. This could, 
of course, happen if the limits for the linear region are set by 
the user of the program to rather small values. The trajectory BE 
illustrates this situation. Should this happen, the trajectory 
will, of course, eventually reach the switching line in the quad- 
rant (j>. pos and negative, where the program would switch 6, again 
back to positive 30 degrees. The system would remain stable and 
eventually again reach the linear control boundary. Of course, 
this would no longer constitute a time-optimal control. By 
choosing the linear boundaries Reasonably large, this situation 
can be avoided. 
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TYPICAL RESPONSES TO COMMAND SEQUENCES 

It was pointed out earlier in this report that an important 
step in the control system design is the examination of the closed 
loop system under a realistic command sequence. Extensive runs 
were made using the program AML75 in Mode I to obtain such response 
time histories. Plotting routines were written to plot the command 
sequence and the associated response. Also plotted, were generally 
the three control surface deflections 6^, 6, and 6 d. 

5 a K 

One typical such response is shown in Figures 18 and 19. The 
run started at level flight at an altitude of 15,000 feet and a 
Mach number of .6. The commanded angles of attack given in one- 
second intervals, reading from left to right, was 

ALPCOM 


8.54 

8.53 

10.70 

8.53 

8.52 

8.52 

12.35 

10.08 

6.56 

5.84 

8.51 

4.74 

16.28 

16.33 

16.63 

8.51 

10.58 

12.19 

13.51 

8.52 

10.75 

12.47 

8.53 

8.51 

8.50 

8.49 

8.49 

10.40 

8.48 

-3.86 

15.94 

while 

8.45 10.15 

the commanded 

11.55 12.7- 8.48 8.48 10.43 

roll angle followed the sequence 

11.89 

-3.07 

PHICOM 

61.15 

61.19 

70.95 

71.72 

72.01 

72.45 

85.48 

91.95 

91.11 

-105 

-117.8 

0 

-64.82 

-74.63 

-84.58 

-57.54 

-62.99 

-65.28 

-65.51 

-49. 

-56.59 

-59.22 

-56.81 

-47.80 

-49.24 

-50.53 

-51.69 

-57.34 

-53.36 

180. 

97.84 

77.34 

80.54 

80.33 

78.63 

6.09 

-165.9 

-166.7 

-167.1 

180. 


Note the transition between nine seconds (4)^ = 91.11) and 10 
seconds ((|)^ = -105). In the plot of the commanded angles, these 
two points are connected directly, the control-system, however, 
executed the command using a right roll, going with 4> through 180 
degrees . 
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- PITCH COMMAND SEQUENCE RESPONSE 
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10.0 20.0 30.0 iO.D 50.0 60.0 70.0 80.0 

TIME, SEC 









Figure 19.- ROLL COMMAND SEQUENCE RESPONSE 
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10.0 20.0 30.0 40.0 50.0 60.0 70.0 60.0 

TIME, SEC 


PART III: DESCRIPTION OF THE 

COMPUTER PROGRAMS 
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PROGRAM AML75 


General Description 

This program represents the main end product of this study. 

It combines the adaptive maneuvering logic with the flight control 
system described in the previous sections of this report. It has 
three basic modes of operation which are summarized below: 

1 ♦ Exercising the program with inputted command sequences. - 
In this mode, the program accepts as input initial conditions for 
the target aircraft and a sequence of commands for roll angles and 
angles of attack. These commands are interpreted by the program 
to be equally spaced in time one second apart. Running the pro- 
gram in this mode is equivalent to exercising the control system 
controlling the nonlinear equations of motion, and, therefore, is 
most suitable for checking out modifications and improvements of 
the control system. 

2. Exercising the program against canned trajectories of the 
attacker .- The purpose of this mode is to have the AML program, 
with its equations of motion being governed by the new control 
system, react against a canned attacker trajectory which is read 
from a magnetic tape. Running in this mode gives the user an op- 
portunity to compare the performance of the AML program with the 
new control system against the old program. This mode also pro- 
vides a convenient tool for debugging the AML program and to check 
out changes in the logic of the tactical decision process. To 
make the program run fast enough so it may be used as a real time 
opponent on the DMS, this mode of running against a canned tra- 
jectory should also prove to be very valuable. 

3 . Exercising the AML program with control system against 
the old AML program .- Here, motion of the target is modeled by the 
same nonlinear equations of motion as on the DMS, while the equa- 
tions of motion for the attacker are those used in the earlier ver- 
sion of the AML program. The purpose of this off-line program is 
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two-fold. First, it demonstrates the capability of the control 
system to "fly" the aircraft effectively under realistic air com- 
bat conditions and to properly execute the AML commanded maneuvers. 
Second, it serves as a stepping stone for later implementation on 
the DMS. In the design of this program, great care was taken to 
retain as much as possible the logic of communication between the 
driver program and the AML program as it is presently in use on 
the DMS. Specifically, the AML program is called in the same man- 
ner as in the present DMS program, that is, with a single argument 
which indicates the function which the AML program has to perform, 
namely, initialization, resetting or (real time) operation. Most 
important, the calling program uses a stepsize comparable to the 
one used on the DMS and calls, during the (real ti me ) execution the 
AML program once every integration step. Internally, however, the 
AML equations of motion use a reduced stepsize. It was found that, 
by using Euler integration, a stepsize of 1/64 of a second was 
sufficiently small to keep the truncation error in the numerical 
solution of the differential equations acceptably small. The 
driver program uses a stepsize of 1/16 of a second (this is the 
stepsize generally used in the AML off-line program), thus, the 
nonlinear equations of motion of the target aircraft and the flight 
control system are integrated over four steps before the AML pro- 
gram returns control to the calling program. 

Structure of the Program AML75 

Figure 20 represents graphically the subroutine structure of 
the AML75 program. Every line from left to right indicates that 
the subroutine to the left of the line has one or several calls 
to the subroutine to the right end of the line. It may appear to 
the casual observer that the structure of this program is unneces- 
sarily complicated. This is, indeed, the case* but it will be 
excused if the history of the development of this program is under- 
stood. 
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The main objective of this contract effort was to design and 
to develop a flight control system accepting as input commanded 
bank angles and commanded angles of attack and then to move the 
aircraft as efficiently as possible from its present attitude to 
the desired attitude. 

The first step in such a design was to select some flight 
conditions which appeared to be fairly representative of air combat 
maneuvering. For a linearized analysis of these flight conditions 
and elementary control loops, stability and control derivatives as 
a function of Mach number, altitude and angle of attack wererre- 
quired. Thus, the first set of subroutines, whose primary purpose 
was to interpolate desired coefficients, were written and debugged. 
Most of this work was carried out on a small, inexpensive computer, 
a Honeywell 716. Program parts from this early phase are recog- 
nizable by their input-output unit assignments (3 and 7 for input, 

4 for output). The subroutines AERFUN and the interpolation rou- 
tines FUNGI, FUNC2A, FUNC2D, FUNC2H, FUNC2X and FUNC3 were devel- 
oped during that first phase. 

The next step consisted of developing the individual programs 
for the linear analysis. These programs are discussed in other 
sections of this report. Of importance here is only the fact that 
ATDYN uses some routines also used by AML75, namely: AERFUN, FUNGI, 

FUHC2A. FUtCfX, and FWMCir 

The lofical next step was to verify that the control system 

based on the linearized equations would also be capable of ad- 
equately controlling the aircraft described by the complete set 

Of nonlinear efeatiens of motion. This was accomplished by ex- 
ercising the AML7S program in Mode 1. The subroutine CMDSEO 

fMntwn ) f r fer o w etsMtiaTly tlHt tttk. ft readt efl 

aerodynamic data and initial conditions for the aircraft (position 
and attitude and their derivatives) in the inertial reference 
frame and a string of angle of attack and bank angle commands and 
then solves the equations of motion for these initial conditions 
and the command sequence. For this program, all the necessary 
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subroutines to integrate the equations of motion were developed, 
such as: 

AERCF: To calculate the coefficients C^, C^, C^, C^, . 

EQMOT: To calculate forces and moments and to calculate u, v, 

w, p, q, f. 

L&GNTR: To control aileron and rudder. 

LOCNTR: To control stabilator. 

DIRCOS: To calculate direction cosine matrix and Euler angles 

from quaternions. 

NLPRNT : To print in condensed form aircraft data. 

Note that so far in the development of the programs, only one 
aircraft was involved. Up to this point, there was therefore no 
need to identify the pertinent aircraft variables as pertaining to 
the attacker or to the target. The subroutines CMDSEQ and the six 
subroutines listed above identify all aircraft variables with 
neither a subscript nor an A or T appended to them. 

After use of the program AML75, running in Mode T, confidence 
in the proper working of the control system was obtained, £he next 
step consisted of providing a command sequence actually generated 
by the AMI. program. The most convenient method to accomplish this 
was to exercise the AML driven program against a canned trajectory. 
This then required interfacing of the AML program and the aircraft 
driven by the control system. Basically, the formerly used sub- 
routine EQMOTT had to be replaced by the sequence of interpolating 
stability and control derivatives, determining rudder, aileron, 
and stabilator deflection, calculating forces and moments acting 
on the aircraft and then calculating u, v, w, p, q and f, and, of 
course, of properly integrating these variables such that position 
and attitude of the aircraft in inertial space was obtained. But 
note, that this new sequence of computations for solving the equa- 
tions of motion had already been accomplished in the previous step. 


78 


except for one thing. When running AML75 in Mode 1, the thrust 
was assumed to be full afterburner thrust at all times, for the 
AML driven aircraft, this is no longer the case. Added to the 
angle of attack and bank angle command was then a thrust command. 
(To keep the program relatively simple, thrust commands of either 
afterburner or idle only were introduced.) To accomplish this step 
three major subroutines were added: RTMAIN3, AMLVS3 and NONLIN. 

The intent was to structure the program so that it could be trans- 
ferred with a minimum of changes to the DMS. The program RTMAIN3 
performs essentially the same function as the program RTMAIN on 
the DMS. Specifically, as RTMAIN on the DMS calls AMLVS2(K) with 
K=1 for initialization, K=2 for resetting and K=3 for real time 
operation, the routine RTMAIN3 calls the subroutine AMLVS3 the 
these three arg.uments for performing the same three functions (the 
operate function is, of course, not a real time function here). 

The trajectory of the attacker aircraft is read from magnetic tape 
every iteration cycle (1/16 of a second) of RTMAIN3. 

In December 1975, the final step of the program development 
was initiated,. Its purpose was to replace the attacker's read-in 
trajectory by a tr of an attacker aircraft driven and com- 

puted by the old equations of motion and the AML program. This 
required the addition of such routines as EOMOTA, THRTLA, REACTA, 
TRYNXA, STATEA and EXTRA. 
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DESCRIPTION OF THE INDIVIDUAL SUBROUTINES 
OF THE AML75 PROGRAM 


Program AML75 (Main Program) 

The program AML75 has three different modes of operation. It 
may simulate the aircraft/control system configuration under a se- 
quence of commands input by data cards (Mode 1). In Mode 2, it 
simulates the AML-driven target aircraft, equipped with the control 
system, maneuvering against a canned trajectory of an attacker. In 
its third mode, it simulates a dogfight between an attacker des- 
cribed by the old equations of motion against a target described by 
the new equations of motion. 

The main program AML75 serves simply as driver for selecting 
the mode. A data card containing the mode is read and then, if it 
is Mode 1, the subroutine CMDSEQ is called, if it is Mode 2 or 3 
the subroutine RTSUB is called. (See Figure 21) 


Subroutine AERCF (Aerodynamic Coefficients) 

Purpose. - This subroutine calculates the nondi mens i onal aero- 
dynamic and moment coefficients. 

Input* - Stability and control derivatives, control surface de- 
flections, weight, C 23 . 

Output* - The three force coefficients C , C , C and the 

X y z 

three moment coefficients C^* C^^, C^^. 

Method. - The following relationships are programmed: 


where Dp, was set to zero and Cp, is obtained by i nterool ati on of 

Uq u 

function F210 (C^^ = f(CL» M)). 
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Figure 21.- Flowchart of Program AML75 
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Subroutine AERFUN 

(Aerodynamic Function Interpolation) 


Purpose ,- To interpolate from the aerodynamic tables the 


s tabi 1 i ty 

and control 

deri vati ves 

and thrust. 


Input 

. - Altitude 

• 3 

Mach number 

, angle of attack. 

Output.- The followinq table 

summarizes the 

functions interpo 

lated by the subroutine 

AERFUN. 



Function 

Symbol 

Fortran 

Name 

Independent 

Variables 

Obtained from 
Tabular Data 
in Array 

Interpol ati on 
Routi ne 
Used 

^L.b 

CLIFTB 

M, 

a» h 

FlOl 

FUNC3 

■-dsB 

CLDLSB 

M, 

a 

Fill 

FUNC2A 

^L. 

a 

CLADOT 

M, 

h 

FI 12 

FUNC2H 


CLQ 

M, 

h 

F113 

FUNC2H 

‘"'-6S 

CLDELS 

M, 

a , h 

F114 

FUNC2H 

6a 

CLDELA 

M, 

h 

F115 

FUNC2H 

^6sp 

GLDLSP 

M 


F116 

FUNCl 


DCDLFT 

M, 

T ri m 

F210 

FUNC2D 

'“asB 

DLCDSB 

M 


F212 

FUNCl 


CYB 

M, 

a, h 

F301 

FUNC3 


CYP 

M, 

a , h 

F311 

FUNC3 


CYR 

M, 

a, h 

F312 

FUNC3 

c 

'6a 

CYDELA 

M 


F313 

FUNCl 

Cy 

6sp 

CYDLSP 

M 


F314 

FUNCl 
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Functi on 
Symbol 

Fortran 

Name 

Independent 
Vari abl es 

Obtained from 
Tabular Data 
in Array 

Interpol ati on 
Routi ne 
Used 

Cy 

''6R 

CYDELR 

M, h 

F315 

FUNC2H 

6 

CLB 

M, a, h 

F401 

FUNC3 

P 

CLP 

M, a, h 

F411 

FUNC3 


CLR 

M, a 

F412 

FUNC2A 

^ t 
^6 a 

CLRDLA 

M , a , h 

F413 

FUNC3 

^6sp 

CLRDSP 

M, a, h 

F414 

FUNC3 


CLRDLR 

M, a , h 

F415 

FUNC3 

'^"'ssb 

CMDLSB 

M, a 

F512 

FUNC2A 


CMQ 

M ‘ 

F513 

FUNCl 

m • 

a 

CMADOT 

M, h 

F514 

FUNC2H 

'tv; 

C|R|LS, . , 

^ * 'ti 


i - J'5 1 5 , 

FUNC3 


CMDELA 

M, a, h 

F516 

FUNC3 

Cm 

•^JCo r. 

6s p 

CMDLSP 

M, a 

F517 

FUNC2A 

'"6 

CNB 

M , a , h 

F601 

FUNC3 

^"P 

CNP 

M , a , h 

F611 

FUNC3 

*^"r 

CNR 

M, a, h 

F612 

FUNC3 


CNDLSP 

M, a 

F613 

FUNC2A 

'"aa 

CNDELA 

M , a j h 

F614 

FUNC3 

“^"aR 

CNDELR 

M, h 

F615 

FUNC2H 
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Fortran 

Name 


Interpol ati on 
Routi ne 
Used 


Functi on 
Symbol 


Obtained from 

Independent Tabular Data 
Variables in Array 


T TMmST H. k 


Fin Of Fill* 


FUMC2H 


Wttkoi.- TtivYar iata ara fivaa far wiif aafla af attack. 


subroutine obtains wing angle of attack by adilaf ana iaarta to 


Tka 


body angle of attack. 


For the coefficient of drag, C. is 

Trim 

C| = C, . + C, 6s + C, l6a| + C. |5sp| 

^Trim ‘‘ds ^6a ^6sp 


Subroutine AMLVS3 (AML Driver for DMS) 

Purpose . - This subroutine will provide the replacement for 
AMLVS2 if and when the new equations of motion are programmed on 
the DMS. It is written in such a manner that the DMS subroutine 
RTMAIN should require no changes. RTMAIN will still call, as in 
the old version, AMLVS. With an argument 1 for initialization, 2 
for resetting and 3 for performing one integration step during the 
operate cycle. 

Input and Output .- The AMLVS3 program during initialization 
reads the aerodynamic tables used by the AML program and during 
operation transfers variables between the subroutine RTSUB (cor- 
responds to RTMAIN on the DMS) and the AML subroutines. 

Method .- A flow chart of the major operations performed by 
AMLVS3 is provided in Figure 22. 


Subroutine CMDSEQ (Command Sequence) 

Purpose. - This is essentially the driver program used to ex- 
ercise the nonlinear equations of motion and the control system if 
a sequence of maneuver commands is read in from data cards. 


♦Table «11 is used for calculation of afterburner thrust. Table 
F831 for idle thrust 
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RETURN 


Kt I UKN 


Fiqure 22.- Flowchart of Subroutine AMLVS3 






OPERATE 

LOOP 



Finure 22 (Continued) 
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To solve the nonlinear 
equations of motion 
for the target with 
a and (j> command and 
thrust command as 



Figure 22. (Concluded) 
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Input and Output .- This subroutine reads all the necessary 
aerodynamic tables for the new equations of motion, the command 
sequence, initial conditions, run parameters and control system 
parameters. It prints trim conditions and, by calling NLPRNT, 
prints the time histories of the flight executing the commanded 
maneuvers. It also writes a magnetic tape containing state vari- 
ables for subsequent plotting by the CALCOMP plotter. 

Method .- Figure 23 shows in form of a flowchart the major 
operations performed by this subroutine. 

Subroutine CSRHOX 
(Speed of Sound, Air Density) 

This is a replicate of the subroutine CSRHO of the AML pro- 
gram, the only difference lies in the storage allocation of the 
tables from which speed of sound and air density are interpolated. 
In the AML program, the two tables for these variables are in the 
common block "TABLES", where as in the new equations of motion 
routines there is a separate common block named "ATMOS". When the 
new equations are programmed to run on the DMS , the DMS functions 
for standard atmosphere should be used rather than reading separate 
tables for the AML program. 

Subroutine DIRCOS 

(Direction Cosines from Quaternions) 

Purpose .- To calculate the direction cosine matrix and the 
Euler angles from the quaternions. 

Input .- The current values of the target quaternions. 

Output. - Target direction cosine matrix, Euler angles. 

Method . - See write-up of subroutine Quatex and Oiler in Refer- 
ence 3. 
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Figure 23.- Flowchart of CMDSEQ 


90 








I 




Figure 23.- (Continued) 


I 
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Figure 23 (Continued) 
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Figure 23 (Concluded) 
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Subroutine EQMOT 
(New Equations of Motion) 

Purpose .- To calculate the angular and linear accelerations. 

Input . - Current values of u, v, w, p, q and r. ^S, the non- 
dimensional coefficients Cy, C^, Cj^, Cj^. 

Method . - Forces are first calculated in the stability axis 
system as 

Xs = ^ S Cx 

Ys = q S Cy 

Z3 = q S C2 

and then transformed to the body axis system 



The moments are obtained directly as 

L ' = q S b 
M = q S c 
N = q S b 

Next, p, q and f, as wel l as u, v and w, are calculated ac- 
cording to the equations derived in the section, "The General Six 
Degrees of Freedom Equations of Motion of a Rigid Body" in this 
report. 
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Subroutines FUNC.. 

(Interpolation Routines) 

Purpose . - To interpolate aerodynamic functions. 

Note ♦ - In the above title, the tv^o periods stand for the dif- 
ferent number and letter combinations which make up the set of 
functions, fUNCT, FUNC2A, FUNC2D, FUNC2H, FUNC2X, FUNC3. The num- 
ber indicates the number of independent variables. The letter, 
which is used in the subroutines with two independent variables, 
indicates what the two independent variables are: 

A Mach, Alpha 

D Mach, Coefficient of lift 

H Mach, Altitude 

X Mach, Altitude for thrust tables. 

Input . - The.jValues of the independent variables to be used in 
the irrterpolatioh.' The function table to be used for interpolation. 

Out p ut . - The interpolat e d function valu e as — a function of 1 ~; 

2 or 3 independent variables. 

Method .- Linear interpolation is used in all routines. If the 
values of the independent variable fall outside the range for which 
the function is tabulated, the closest value in the table is used 
(no extrapolation). In FUNC2H, interpolation for altitude is only 
performed if h lies between 15,000 feet and 45,000 feet. For 
h £ 15,000 feet, the function value at 15,000 feet is taken, for 
h ^ 45,000 feet, the function value at 45,000 feet is used. 

Subroutine GETCOM (Get Command) 

Purpose . - To translate the AML commands into commands for bank 
angle and angle of attack. 

Input .- "New command" flag (ICMNWT), AML commanded load factor 
(as fraction of g^^^ available), AML command maneuverplane rotation 
angle CROTT). 
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Output .- Commanded bank angle (PHICMD), commanded angle of 
attack (ALFCMD). 

Method . - First, the subroutine checks whether AML has command- 
ed a new maneuver; if not, the old and remain unchanged. 

If AML has issued a new command, the maximum available loadfactor 
(FLODMT) Is calculated (by a call to AERF4). The actual loadfac- 
tor is then the product of desired load factor level and the maxi- 
mum load factor available. The angle of attack may be calculated 
as follows: 


L = w ♦ load factor 

a = C, 

qS 

The. commanded bank angle is obtained from the AHL c^^^ 
maneu verpl ane rotation angle by noting that the aircraft must be 
rotated by an angle c|)* about the maneuverpl ane x-axis. This will 
create a component of the lift to compensatfiv for the gravity force 
along the maneuverpl ane y-axis. The procedure for calculating (()* is 
explained iri detail on page 19 and following of . Reference 2. 

Subroutine LACNTR (Lateral Control) 

Purpose . - To determine aileron and rudder deflections. 

Input .- Aircraft state variables, control system parameters. 

Output .- Aileron deflection (6a) and rudder deflection (6R). 

Method . - The functional relationships used are described in 
the earlier section, "The Lateral Control Systems". 
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Subroutine LOCNTR 
(Longitudinal Control) 

Purpose .- To determine stabilator deflection to achieve com- 
manded angle of attack. 

Input .- Aircraft state variables, control system parameters, 
commanded angle of attack. 

Output . - Stabilator deflection 6^. 

Method . - Basically, the control system as described in the 
section "The Longitudinal Control System" is implemented. 

If Kp is zero, no longitudinal control is performed. Other- 
wise, a and ^ are computed first. 

The following gains for q < 150 Ibs/sq. inch are set: 

Kp = 3 Kj = 3 Kp = 0.6 

a a a 

For dynamic pressure q > 520 Ibs/sq. inch, the following gains 
are scheduled: 

K = 2 Kj = 4 Kp = 0.25 

For dynamic pressures between 150 Ibs/sq. inch and 520 Ibs/sq. 
inch, the gains as read in from the input card are used. 

Next, the commanded angle of attack is reduced if the present 
roll angle error is greater than 30 degrees and the roll rate ex- 
ceeds 50 degrees/second. 

Stabilator deflection according to Figure 12 is next calcu- 
1 ated . 

Finally, stabilator deflection is limited to -23 degrees and 
+7 degrees. 
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Subroutine NONLIN 
(Nonlinear Equations of Motion) 


Purpose . - This subroutine is the driver for the routines per- 
forming the solution of the nonlinear equations of motion in Mode 2 
and Mode 3 of AML75. 

Input and Output .- The aerodynamic data for the target air- 
craft are read by this subroutine during initialization. During 
the operate mode, the subroutine effects printout of some aircraft 
key variables by calling NLPRNT, it also transfers target data 
through the common block OPDATl to the calling program (AMLVS3). 

It receives maneuver commands from AMLVS3 through common block 
0PDAT2. 

Method . - It is important to realize that the only communica- 
tion between NONLIN and AMLVS3 is through the two common blocks 
OPDATl and 0PDAT2. Also note that NONLIN uses a different inte- 
gration stepsize than the AML program. However, in both programs 
the Fortran name for the integration step is DT ; the mechanism, by 
which these two different variables are kept separate is by having 
one of them in common block K4, the other one in common block K5. 
All the routines associated with the AML program contain common K4, 
all routines associated with the solution of the nonlinear equa- 
tions of motion contain common K5. 

Otherwise, the subroutine NONLIN is very similar to the sub- 
routine CMDSEQ, of which a flowchart is provided in Figure 
The main differences between NONLIN and CMDSEQ are: NONLIN is 

called with one argument, which is 1 for initialization, 2 for re- 
setting, and 3 for operation, whereas in CMDSEQ these three func- 
tions are performed in sequence within CMDSEQ. NONLIN does not 
read the command sequence data deck. NONLIN transfers aircraft 
state variables to common OPDATl. NONLIN picks up commands from 
common 0PDAT2. 
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Subroutine NLPRNT 
(Print for Nonlinear Program) 

This subroutine prints the following target variables when- 
ever it is called by NONLIN or by CMDSEQ. 

^ ^ ^ ^ '^Tot ^ 6a 6R (SS e 4 ) h 

All angles are printed in degrees, angular rates in degrees/second. 

Subroutine PAIRCR 
(Print Aircraft Data) 

This subroutine corresponds to the subroutine PRNTAC of Ref- 
erence 3 and produces printout as illustrated in Figure 11 of 
Reference 3. At the present time, not all variables printed by 
this subroutine may have been transferred to the corresponding 
common block from which PAIRCR selects the variables to print. 
(Example: P, Q and R for target, but these variables are printed 
in NLPRNT. The time interval between calls to PAIRCR is set to a 
fixed value of one second in subroutine RTSUB. 

Subroutine PCELL 
(Print Tactical (Cell) Data) 

The printout generated by this routine allows for quick assess- 
ment of the relative tactical situation between the two aircraft. 

The printout is very similar to the one shown in Figure 12 of Ref- 
erence 3, the annotations to questions 5 through 13 do not corres- 
pond to the actual questions analyzed in the routines STATEA and 
STATET. 


Subroutine RTSUB 
(Real Time Subroutine) 

This subroutine simulates the functions performed on the DMS by 
the subroutine RTMAIN. For details on input, output and method, see 
flowchart of RTSUB, Figure 24. 
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Figure 24 (Continued) 
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OPERATE 

LOOP 


Figure 24 (Concluded) 
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Trajectory 


GO TO 
IMODE 


Against AML 


READ 1 RECORD OF 
ATTACKER DATA 
FROM MAG TAPE 


CALL 

EQMOTA 


CALCULATE 
ATTACKER 
DIRECTION 
COSINE MATRIX 


CALCULATE 
ATTACKER 
DIRECTION 
COSINE MATRIX 


CALL 

AMLVS3(3) 


CALL 

AMLVS3(3) 


T IME 

r^^TO PRINT?. 


CALL 

REACTA 


CALL 

STATET, PRCELT, 
PAIRCR, PCELL 


CALL 

THRTLA 
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TO PRINTS 


^MORE \ 
POINTS ON 
\ TAPE?^ 


CALL RELGN, 
;TATEA, STATET, 
•RCELT, PAIRCR, 
PCELL 


STOP 


NO YES 

POINTS TO 
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Figure 24.- (Concluded) 














INPUT DATA FOR AML75 


The input data for this program depend on the mode in Which 
the program is used. For evaluation of command sequences (Mode 1) 
no input data associated with the AML program is required. Figure 
25 gives an overview of the input data deck for Mode 1. In Modes 
2 and 3, data for both the AML program and the new equations of 
motion are required. The input data decks for Modes 2 and 3 are 
the same and are shown in Figure 26. 

Certain blocks of data are used in different places and will 
therefore be described only once. 

Input Data for Mode 1 Operation 

Program Mode .- First card of deck specifies program mode. 

This must be a 1 for Mode 1 operation. 

1 Card FORMAT (110) 

Aerodynamic Functions for New Equations of Motion .- In all 
the folldiMng tables, unless specified otherwise, data listed as 
function of Mach number are given for the following 7 Mach numbers: 

0.2, 0.4, 0.6, 0.8, 0.9, 1.0, and 1.1. 

Functions tabulated as function of angle of attack are given 
for the following angle of attack values: 

-5°, 0°, 5°, 10°, 15°, 20°, 25°, and 30°. 

Functions tabulated as function of altitude are given for the two 
altitudes 15,000 feet and 45, 000. feet. 

(a) Basic coefficient of lift ( C^ basic’ Fi**iction FlOl) is tabu- 
lated as function of Mach, angle of attack and altitude. 

16 Cards FORMAT (7F10.7) 
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INPUT DATA TO AML75 
Mode 1 (Command Sequence) 



Figure 25 . -Input Deck Setup for AML75 (Mode 1) 




Read in 
RTSUB 


Read in 
NONLIN (1 


Read in 
AMLVS3(1 ) 


Read in 
AML75 
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I 


(b) Coefficient of lift with respect to speedbrake deflection 

(C. , Function Fill) is listed as a function of Mach and al- 

6sB 

ti tude . 


8 cards FORMAT (7F10.7) 

! (c) Coefficient of lift with respect to g Function F112) is 

listed as function of Mach and altitude. “ 

2 cards FORMAT (7F10.7) 

I 

(d) Coefficient of lift with respect to q , Function F113) is 
listed as function of Mach and altitude. ^ 

I 

2 cards FORMAT (7F10.7) 

I ^ Coefficient of lift with respect to stabilator deflection 

(C| , Function F114) is listed as function of Mach, angle of 

6s 

attack and altitude. 


16 cards FORMAT (7F10.7) 

(f) Coefficient of lift with respect to aileron deflection (C|^ , 

F115) is listed as function of Mach and altitude. 

2 cards FORMAT (7F10.7) 

(g) Coefficient of lift with respect to spoiler deflection (C^ , 

F116) is listed as function of Mach 

1 card FORMAT (7F10.7) 

(h) Coefficient of drag due to lift (C^, F210) is listed as func- 

tion of Mach, and coefficient of lift for the following 13 
values of C|_: 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 1.0, 

1.2, 1.4, 1.6 


(i) 


13 cards 


FORMAT (7F10.7) 


Coefficient of drag due to speed brake deflection 
is listed as function of Mach. 


(Cn . F212) 
■^sb 


1 card 


FORMAT (7 FI 0*7) 
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Side-force derivative with respect to side slip angle (Cy , 

F301 ) is listed as function of Mach, angle of attack and alti- 
tude . 

16 cards FORMAT (7F10.7) 

Side-force derivative with respect to roll rate (C^ , F311) is 

listed as function Of Mach, angle of attack; altitude. 

16 cards FORMAT (7F10.7) 

Side-force derivative with respect to yaw rate (Cw , F312) is 

^ r 

listed as function of Mach, angle of attack and altitude. 

16 cards FORMAT (7F10.7) 

Side-force derivative with respect to aileron deflection 

(Cw , F313) is listed as function of Mach. 

6a 

1 card FORMAT (7F10.7) 

Side-force derivative with respect to spoiler deflection 

(Cw 9 F315) is listed as function of Mach and altitude. 

6sp 

1 card FORMAT (7F10.7) 

Side-force derivative with respect to rudder deflection 

(Cw , F315) is listed as function of Mach and altitude. 

^6R 

2 cards FORMAT (7F10.7) 

RolVijlg-iomeint derivative with respect to side-slip (G« , 

F401 ) is listed as function of Mach, angle of attack and al- 
ti tude . 

T6 cards FORMAT (7F10.7) 

Rolling-moment derivative with resoect to roll rate fC» . 

F411) is listed as function of Mach, angle of attack and al- 
t1 tuMie . ' ■ . , ’ - 

FORMAT (7F10.7) 


16 cards 


Rolling-moment derivative with respect to yaw rate (C„ , F412) 
— - - _ 

is listed as function of Mach and angle of attack. 

8 cards FORMAT (7F10.7) 

Rolling-moment derivative with respect to aileron deflection 

(Cp , F413) is listed as function of Mach, angle of attack 
^6a 

and altitude. 

16 cards FORMAT (7F10.7) 

Rolling-moment derivative with respect to spoiler deflection 

(Cp , F414) is listed as function of Mach, angle of attack 
^6sp 

and altitude. 

16 cards FORMAT (7F10.7) 

Rolling-moment derivative with respect to rudder deflection 

(Cp , F415) is listed as function of Mach, angle of attack 
■^6R 

and altitude. 

16 cards FORMAT (7F10.7) 

Pitching-moment coefficient (C^ 1 ^, F501 ) is listed as func- 
tion of Mach, angle of attack and altitude 

16 cards FORMAT (7F10.7) 

Pitching-moment derivative with respect to speed brake de - 
flect! on (C , F512) is listed as function of Mach and angle 
of attack. 

8 cards FORMAT (7F10.7) 

Pitching-moment derivative with respect to pitching rate {0^^ , 
F513) is listed as function of Mach number. 

1 card FORMAT (7F10.7) 

Pitching-moment derivative with respect to a F514) is 

listed as function of Mach and altitude. “ 

FORMAT (7F10.7) 


2 cards 




8 cards FORMAT (7F10.7] 
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(gg) Yawing-moment derivative with respect to aileron deflection 

(C , F614) is listed as function of Mach, angle of attack 
6a 

and altitude. 

16 cards FORMAT (7F10.7) 

{ hh ) Yawing-moment derivative with respect to rudder deflection 

(C , F614) is listed as function of Mach and altitude. 

"6R 

2 cards FORMAT (7F10.7) 

(ii) Afterburner thrust (F811) is listed as function of Mach and 
altitude for the following altitudes (sea level# 10k, 20k, 
30k, 40k, 50k, 60k feet) 

7 cards FORMAT (7F10.7) 

(jj) Idle thrust (F831) is listed as function of Mach and altitude 
for the same values of altitude as in (ii) above 

7 cards FORMAT (7F10.7) 

The entire data deck for the aerodynamic and propulsion data 

for the new equations of motion is composed of 356 cards . 

Atmospheric data .- The next block of data is read into the 
tables in common block ATMOS. It consists of 

(a) Speed of sound in feet/sec is tabulated as function of alti- 
tude from sea level to 60,000 feet at 500 feet increments. 

12 cards FORMAT (10F8.7 

O 

(b) Air density in slugs/feet is tabulated for the same altitude 
points as (a) above. 

13 cards FORMAT (10F8.7 

The entire atmospheric data block contains 26 cards. 

Command sequence .- The next data block contains the commanded 
values of bank angle and angle of attack. Values are listed in 

pairs of (<t>com* ®com^* ^ ® expressed in degrees , Jthe cQinipands 

are interpreted by the program as being one second apart. The 


no 


first pair of commands is executed at time = 0, the second pair at 
time = 1 second, and so forth. The program expects command input 
for 199 seconds, however, only command input from time zero until 
Tend actually used. The input cards are prepared as follows: 



Time=0 Time=l sec Time=2 sec Time=3 sec 

50 cards FORMAT (8F10.2) 

Title card . - th« card contairiSi a Jcun ti t1§, vi|Ji.iQh is 
printed at the beginning of the printout of the program. 

1 card FORMAT (72A1) 

Initial conditions .- The next card contains the information 
required to establish the initial conditions for the run. The 
following variables are listed on this card: 

x„ inertial x-coordinate 

e 

inertial y-coordinate 
h inertial altitude 

Vfot Magnitude of initial velocity vector 

body axis yaw Euler angle 
0 body axis pitch Euler angle 

4 > body axis roll Euler angle 

a initial angle of attack 

$ initial sideslip angle 

Note that the initial angle of attack and the initial side- 
slip angle may >e. changed when the run requests certain trim op- 
tions. (See section "Use of Options for CMDSEQ".) 

1 card FORMAT (7F10.2, 2F5.2) 

Initial control surface deflections .- Initial aileron, rudder 

and stabilator deflection are listed pn tpis ;cPr4. flote aga4|f 
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that, depending on trim options, these values may be changed by 
the program. The initial spoiler deflection is set by the program 
to 1.4333 — times the initial aileron deflection. Input are in 
degrees . 

1 card FORMAT (3F10.2) 


Run parameters .- This card contains the following variables: 


DT 


T 


end 


KPRINT 


IVARd ) 


IVAR(3) 


Integration stepsize. 

This is the integration stepsize for the solution of 
the nonl i near equations„of motion and it is usually set 
to 1/64" = 0.015625. 

End time. 

Generally, command sequences should be provided for time 
equals zero up to T^^^ minus one second. 

A variable controlling the time interval between calls to 
subroutine NLPRNT. This time interval is equal to 
KPRINT times DT. 

Angle of attack trim control variable. 

If IVAR(l) is set to 1, the initial angle a is set to 
the trim angle of attack for trimmed, level flight. If 
the input initial conditions are such that trimmed level 
flight is not possible, the program will print this fact 
and will set the angle of attack to 30.9 degrees. 

Stabilator trim control variable. 

If this variable is set to one, the program will set 
the initial stabilator deflection to the trim value and 
thus override the previously inputted initial stabila- 
tor deflection. 


IVAR(4) If this variable is set to 1, the program will keep the 
thrust during the entire run equal to the initial trim- 
med value. 

1 dard FORMAT (2Fi 0.6, 4 
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Longitudinal gains .- This card reads the following variables 
and parameters: 

VAR(2) Commanded angle of attack. 

This is a remnant of an earlier version of this 
program (before entire command sequences were read 
in) and is no longer used. The input in this field 
is ignored by the program. 

VAR(3) Longitudinal gain Kp . 

VAR(6) LongitudinaT'gain K. . ’ ' ' . 

■ " ' " ''' ’ ' ' ^ 

VAR(IO) Longitudinal gain Kp (a- rate gain) 

a 

These three gains are used by the program during flight con- 
ditions of low dynamic pressure £ 840 psi). The high dynamic 
pressure gains are determined in the subroutine LOCNTR. 

1 card FORMAT (4F10.3) 

Lateral gains and control parameters .- The following param- 
eters associated with the lateral control system are listed on this 
card: 

VAR(4) Commanded bank angle. 

The same comment as for VAR(2) in the preceding 
card applies here. 

VAR(5) Roll position gain Kp . 

4 > 

VAR(7) Roll integral gain Kt . 

4 > 

VAR(9) Roll rate gain Kp . 

VAR(8) Beta position gain Kp . 

8 

6N(1) Beta rate gain Kn . . 
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GN(2) Bang-bang position threshold. 

GN(3) Bang-bang rate threshold. 

The gains listed here are used for low q flight conditions 
(q < 840 Ibs/square inch), gains for high dynamic pressure are de- 
termined in subroutine LACNTR. 

1 card FORMAT (8F10.3) 

Use of Options for Mode 1 

By cKSolift^^ c of input options' a number of differ- 

ent types of controlled or uncontrolled responses can be obtained 
from CMDSEQ. 

The flight condition is established by altitude, absolute 
(HAB) and total velocity (VTOTAL). With the body-axis Euler angles 
, 6, and (j) any aircraft attitude initial condition can be input. 
Initial angle of attack a (ALPDEG) and initial sideslip ange 
(BETZER) establish the orientation of the relative wind with 
respect to the body axes. 

The second card sets the initial deflection of aileron 6^ 
(DLADEG), rudder 6^^ (DLRDEG), and stabilator 6^ (DLSDEG). Refer- 
ring to card 4, if VAR(3) = Kp is set to zero the stabilator de- 

6 

flection 6^ will be left at the initial value during the run. 

Likewise, if VAR(5) = Kp is set to zero, the aileron deflection 6 

<|5 ^ 

and the rudder deflection 4^ will both be leftat their initial 

settings. This provides a convenient way of making the aircraft 

fly at fixed control surface deflections, be they zero or nonzero. 

The trim control integers of card 3 cause the program to cal- 
culate and set the corresponding variable initially at the trimmed 

level flight value. Thus, if IVAR(l) = is set to 1 , the 

a 

program will calculate the trim angle of attack and set a(ALPDEG) 

initially at this value. Likewise. IVAR(3) = 1^1^ =lwill , 

stab 

result in setting the stabilator initially at the trim value and 
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IVAB(4) - tiirust level at the trimmed 

thrust 

value. The only way to keep the stabilator at the trimmed value is 
to set VAR(3) = Kp =0. 

a 

For an example, the following settings would result in uncon- 
trolled flight starting at a trimmed condition: 

IVAR(l) = IVAR(3) = IVAR(4) = 1 

VAR(3) = VAR(5) = 0 

Input values of ALPDEG from card 1 and DLSDEG from card 2 would 
have no effect. 

However, the following settings would result in uncontrolled 
but initially trimmed flight for the longitudinal axis and control- 
led flight in the lateral axes. 

IVAR(l) = IVAR(3) = IVAR(4) = 1 

VAR(3) = 0, VAR(5) H 0. 

The desired gains for roll control would be input via card 5. 

Input Data for Modes 2 and 3 Operation 

A complete data deck setup for Mode 3 is shown in Figure 26 
and Appendix B lists a complete set of input data required for 
Mode 3. 

Program mode .- The first card contains the specification for 
the desired program mode. Mode 2 requests exercising the AML pro- 
gram against a canned trajectory. Mode 3 against an AML program 
using the old equations of motion. 

1 card FORMAT (110) 

Atmospheric model definition for old equations of motion . - 
The speed of sound and air density tables as explained in the cor- 
responding section of the Mode 1 input data description. 

26 cards FORMAT (10F8.7) 
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Mach number and altitude definitions .? Xn the following tables 
listing functions used for the old equations of the AML program, 
data given as function of Mach number are listed for the following 
10 values: 0.2, 0.5, 0.8, 0.9, 1.0, 1.1, 1.2, 1.5, 1.8 and 2.0. 

Functions listed as function of altitude are given for the follow- 
ing 5 altitudes: sea level, 15k, 30k, 45k and 55k feet. 

Maximum permissible load factor .- The maximum permissible load 
factor in g's is listed as function of Mach and altitude. 

5 cards FORMAT (10F8.7) 

Mi niiBiW ftl lowa#le Mach^* number . - This vari able Is listed as 
function of altitude. 

1 card FORMAT (5F7.3) 

Maximum allowable Mach number .- This variable is listed as 
function of altitude. 

1 card FORMAT (5F7.3) 

Lift coif Hcient derivative with respect to angle of attack 

(C|^ ).- This variable is listed as function of Mach number and 

ot 

al ti tude . 

5 cards FORMAT (10F7.3) 

Military thrust as function of Mach number and altitude . - 
Thrust for one engine is listed as function of Mach number and al- 
titude for the following 14 values of the Mach number: 0.2, 0.4, 
0.6, 0.8, 0.9, 1.0, 1.1, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2 and 2.4, 
and for the following 7 values of altitude: Sea level, 10k, 20k, 

30k, 40k, 50k, 60k feet. 

14 cards FORMAT (7F10.0) 

Afterburner thrust .- Afterburner thrust is listed as function 
of Mach and altitude for the same Mach numbers and altitudes as 
military thrust. 

14 cards FORMAT (7F10.0) 
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Idle thrust *- Idle thrust is Visted as of Mach, and 

altitude for the same data points as military thrust. 

14 cards FORMAT (7F10.0) 

Coefficient of drag .- The coefficient of drag is listed as 

function of coefficient of lift and of Mach number. Drag coeffi- 
cients are given for the following 16 values of C^: 0, 0.1, 0.2, 

0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4 and 1.5. 

30 cards FORMAT (6F8.4/6F8.4/4F8.4) 

Pi ve. recovery angle .- Maximum dive angles, in degrees, from 
which the aircraft can pull out without hitting the ground, are 
listed as functions of Mach number and altitude for the following 
5 values of altitude: Sea level, 3k, 6k, 9k and 12k feet. 

5 cards FORMAT (10F7.3) 

Sustained load factor .- The load factor at which the air- 
craft can perform a turn when in afterburner without losing energy 
is listed as function of Mach and altitude. 

5 cards FORMAT (10F7.3) 

The entire aerodynamic and propulsion deck for the old equations 
of motion consists of 94 cards . 

Aerodynamic tables for target .- These are the tables FlOl 
through F831 as explained in detail in the section "Input Data for 
Mode 1 Operation". 

’ 356 cards FORMAT (7F10.7) 

Atmospheric data .- This data check is identical to the one on 
atmospheric data documented above. 

26 cards FORMAT (10F8.7) 

Longitudinal gains and end time .- This card contains the 
following information: 

VAR(2) Commanded angle of attack (not used 

in this program. 
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VAR(3) 

Proportional a-feedback gain K 

VAR(6) 

Integral a-feedback gain Kj . 

a 

VAR(IO) 

Alphadot feedback gain Kj^ . 

a 

TEND 

End time. 


Notes: 1) The gains are valid between 150 Ibs/sq. inch 530 

Ibs/sq. inch gains below and above these limits are 
de^frmined in subroutine LOCNTR. 

t) There are two values of T^^^ read by this program (the 
one read here should be deleted from the program). The 
smaller value of T^^^ will always terminate the run. 

1 card FORMAT (8F10.3) 

Lateral gains and control parameters .- This card is identical 
to the corresponding card in Mode 1, see documentation "Input data 
for Mode 1 Operation". 

1 card FORMAT (8F10.3) 

End time .- This card contains the end time Tend in seconds. 

1 card FORMAT (lOX, F10.2) 

Target aircraft initial conditions .- The target aircraft in 

Modes 2 and 3 is the one modeled by the new equations of motion. 
This card lists the initial conditions as follows: 

Xg Target initial x coordinate in inertial coordinate 

system . 

y^ Target initial y coordinate in inertial coordinate 

system . 

h Target initial altitude in inertial coordinate sys- 

tem. 

Xg Target initial Xg velocity component. 
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Yg ' T«pg«t liitftlal vel^oclty comp-onent 

fi Target initial! altitude dot component 

4» Target initial body Euler yaw angle in degrees. 

0 Target initial body Euler pitch angle in degrees. 

(p Target initial body Euler roll angle in degrees. 

Note: The subroutine NONLIN resets the target aircraft to trimmed 

level flight so that the input of 0 and (fi are actually ignored. It 
is also . .adyiiablei.iia input B as zero ^VQld conf] 1ct,,wi th the . 
trim cbif^ffdih 

1 card FORMAT (8F10.2) 

Attacker aircraft initial conditions. - This card is required 
only in Mode 3, since in Mode 2, the attacker initial condition is 
read from the magnetic tape containing attacker data. The list of 
variables on this card is the same as the one on the preceding 
card for the target. 

1 card FORMAT (8F10.2) 

Summarized, the entire input data deck for Mode 2 contains 507 
cards, for Mode 3 it contains 508 cards. 
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PROGRAM ATDYN 


General Description 

The program ATDYN was originally envisioned as a self-contained 
program which would evaluate and print out complete information on 
the attitude dynamics of the aircraft for a specified flight con- 
dition. Stability derivatives calculated from tabular data were 
used to find the coefficients of the linearized longitudinal and 
lateral equations of motion. These equations were then to be used 
to find the transfer functions needed in the respective control 
loops of the longitudinal and lateral axes. It was also deemed de- 
sirable to find the coefficients for the first-order differential 
equation form of these equations for use in state variable analysis 
and synthesis techniques. 

All closed-loop linear control analysis capability was incor- 
porated into other programs than ATDYN. However, since a very 
limited pure proportional gain closed-loop root option exists in 
ATDYN, it is documented here for the sake of completeness. 

The first approach used for attitude dynamics was to solve 
algebraically for the coefficients of the polynomial for the de- 
sired characteristic equation and use the subroutine POLFAC to 
factor the resulting polynomial. This approach is exemplified in 
the subroutine LATPOL which calculates the roots of the lateral 
characteristic equation. This approach is direct, reliable, and 
computationally fast but requires tedious algebraic manipulation 
for each polynomial needed. This algebraic approach was used 
again in the subroutine ALFZER to find the numerator roots (zeros) 
of the transfer functions a(s)/6^(s) and a(s)/6^(s) for aileron 
and rudder input, respectively. To find all the transfer functions 
needed using this approach would have been extremely tedious. 

Also, problems with the subroutine POLFAC indicated a need for an 
alternative root-finding routine. Furthermore, a procedure easily 
adaptable to both open and closed-loop analysis was needed. 
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After some experimentation with eigeTivalue routines the -rele- 
vant subroutines from the EISPACK, the eigenvalue subroutine pack- 
age developed by Argonne National Laboratory were applied. These 
routines provided a trustworthy eigenvalue procedure but finding 
numerator dynamics of transfer functions was not convenient. 

Shortly thereafter the subroutine package TRANS was made opera- 
tional. TRANS was developed specifically for control analysis and 
was therefore more convenient than any previously tried technique. 
Thereafter, all dynamics analysis was done using TRANS, but the 
previously «devel oped routines were kept for comparison. 

structure of Program ATDYN 

The program ATDYN consists of a main program and 22 individual 
subroutines. Of these, the subroutine CMTRX is documented in Ref- 
erence 3, the following eight subroutines are documented in the 
section "Description of the Individual Subroutines of AML75" in 
this report: CSRHOX, AERFUN, FUNGI, FUNC2A, FUNC2D, FUNC2H, 

FUNC2X and FUNC3. Figure 27 illustrates the calling hierarchy of 
the program ATDYN. The function of the remaining subroutines is 
briefly explained below. 

Subroutine KS .- This subroutine computes the 19 coefficients 
K| through K-jg used for the linearized longitudinal equations and 
the 22 coefficients K^q through Kg-j used for the lateral linear- 
ized equations. It also punches on cards the entire array of the 
coefficients through . (This deck is used as input to a 
number of the linearized analysis programs.) 

Subroutine ALFZER .- This is a specialized routine to calcu- 
late the numerator polynomial coefficients for the following 
transfer functions 

a(s)/6g(s) (|)(s)/6g(s) d>(s)/6j^(s ) . 

It then finds the zeros of these pdlynomials. 
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BALANC 

ELMHES 

HQR 


POLRT 


Figure 27.- Structure of Program ATDYN 


Input to ALFZER consists of the K-array and a flag NCASE with 
the following function. 

NCASE = 0 Calculates zeros of a(s)/6^(s) 

NCASE = 1 Calculates zeros of (b(s)/6 (s) 

NCASEr 2 Calculates zeros of <j> ( s )/6|^( s ) 

These zeros are printed out in ALFZER and are not returned to 
the calling program. 

Subroutine LATPOL .- This is a specialized subroutine which 
calculates the polynomial coefficients of the lineaHzed lateral 
characteristic equation, and finds the corresponding roots. Input 
to LATPOL is the array of coefficients K^. , the roots are printed 
out by LATPOL but not returned to the calling program. 

Subroutine POLFAC .- This is a general purpose polynomial fac- 
toring routine which sometimes fails to give the correct roots. 

Subroutine FACTOR, POLDRM, POLRT .- These three routines per- 
form the polynomial factorization for this subroutine TRANS. 

Subroutines RG, BALANC, ELHES, HQR .- These are adaptations of 
subroutines from the Eigervalue package EISPACK, developed by 
Argonne National Labs and they are documented in detail in the 
writeup of EISPACK. 

Subroutine TRANS .- This is a general-purpose subroutine for 
calculating transfer functions of linear systems directly from the 
Laplace transformed equations of motion of the system. The use of 
TRANS will be illustrated using the equations for the lateral axis 
as exemplified in LABEZR. 

If the Laplace transforms of the lateral axis linearized 
equations are written in terms of the variables <j)(s), i|^(s), 3(s) 
with forcing functions 6,('s) and 6„(s) we obtain the following 

a I 

matri x : 
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<|)(s) 


3(s) 

4,(S) 

5^ts) 

(KjoS^+KjiS) 

(K32S2+K33S) 

•^34 

*^35 

*^36 

( K^yS^+K^gS) 

(K 35 S^K 4 oS) 

>^41 

K42 

^43 


(K46S+K47) 

(K48S+K49) 

•^50 

^51 


If these linear equations are solved accerding to Cramer's 
rule for tlie variable tpt I) for rudder input 5^(s) we obtain 


Do(s) 6 (s) 

3($) = D(s) where 

















nal equations. This is done using the three-dimensional array G*. 
Thus, C(i, j, k) represents the element of the i-th row and j-th 
column which is the coefficient of s . Starting with the zeroth 
powers of s, we find; 


C(1.3;l) = K 34 

C(l,l,2) 

= 

•^31 

C(l,l,3) = 

So 

C(l,4,l) = K 3 g 

C(l,2,2) 

== 

•^33 

C(l,2,3) = 

S 2 

C(2,3,l) = 

C(2,l,2) 

= 

^ 38 

C(2,l,3) = 

•^37 

e(2.C.tl - K 43 

C(2.,2,2) 

a 

^40 

C(2,2,3) =, 

.K 3 , 

C(3,l vl) -.^5 

C(3,l,2) 

a 

^44 



C{3.2.U = 447 

C(3,2,2) 

= 

^46' 



c(^.3.y >r K 49 

-i c. V .j, ■ 

C(3.3,2) 

= 

48 



- ■■■ ■>; Iw ■ * 

C(3,4,l) = K 5 I 






Since TRANS will 

use only 

one 

1 right 

-hand side. 

if we want 6 ^(s) 

the forcing function we must 

put the 

coef f i ci ents 

of 6 (s) in 


column 4 of the C matrix. 

It would have been better if we could use integer zero to 
represent s° but Fortran does not permit zero to be used as the 
index of an array. 

If the C array is correctly input the subroutine TRANS will 
calculate properly the polynomial D(s). 

To determine the numerator polynomial Dg(s) we wish to replace 
column 3 by column 4 ( 6 ^(s) coefficients) and evaluate the corres- 
ponding determinant. This is done by means of the 5-element array 
NMRTR. To obtain the numerator polynomial corresponding to the 
third column we would let NMRTR = (3, 0,0, 0,0). 

The call to TRANS would then be: 

CALL TRANS(3,2,NMRTR) 
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where the first nuf^er, 3* 1 s the .s,i zei ot or? g matrix, the 
second number, 2, is the maximum power of s in the matrix, and 
NMRTR was described above. 

Input Data Deck Setup for Program ATDYN 

Figure 27 illustrates an input data deck to the program ATDYN. 
The first two blocks of the data deck, consisting of the aerodynam- 
ic functions and the atmosphere data, are explained darlier in this 
report., for each individual case to be analyzed by ATDYN, a "RUN 
DATA” dei||| 5 ^ df five cards is required. These five cards 

contain the following information: . 

- '*■ rA--,.. 

Car^itf T . - fi tl e card. The first 78 ciharacters o| Jhis card ^ 
will be printed at the beginning of the printout. If the title 
card contains END in columns 1 through 3, the computer run will be 
terminated. 

Card 2. - Run options, FORMAT (41.2). These option parameters 

are : 

NOPT: Option number, may be an integer between 0 and 5, spec- 

ifying the following options: 

0 = longitudinal open loop roots only. 

1 = longitudinal roots plus closed loop root with gain 

Kp , specified on card 4. 
a 

2 = lateral open loop roots only. 

3 = lateral open loop roots plus lateral closed loop root 

with gain Kp , specified on card 5. 
a 

4 = open loop longitudinal and lateral roots. 

5 = all options (0 through 4) combined. 

IPTRAN: A parameter controlling polynomial printout in subroutine 

TRANS. 

IPTRAN = 0: Part of polynomial printout will be suppres- 

sed. 
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t 




Figure 27- Input Deck Setup for Program ATDYN 
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IPTRAN = 1: All polynomials will be printed. 

JPTRAN: A parameter controlling root printout in subroutine TRANS. 

JPTRAN = 0: Part of root printout will be suppressed. 

JPTRAN = 1: All roots will be printed. 

KPTRAN: A parameter controlling matrix printout in subroutine 

TRANS. 

KPTRAN = 0: printout of matrix elements will be suppres- 

sed. 

KPTRAN = 1: matrix elements will be printed. 

The most frequently used options when running ATDYN for rou- 
tine control analysis is: 

N0PT=4 IPTRAN=0 JPTRAN=1 KPTRAN=0 

Card 3. - Flight condition, FORMAT (8F10.2). The following 
variables Ihre spectfied oii this card. 

e Inertial x coordinate (not used by program) 

y^ Inertial y coordinate (not used by program) 

h Altitude 

Vfot Magnitude of velocity vector 

ijj Euler yaw angle (not used by program) 

0 Euler pitch angle in degrees 

(j) Euler roll angle in degrees 

a Body angle of attack in degrees 

Card 4 .- a-feedback gains. FORMAT (8F10.2). This card may be 

left blank unless option 1 or option 5 is used. In this case, 

eight values of gains Kp must be specified. The program will cal- 

a 

culate the closed loop poles for these eight gains. However, since 
ATDYN has no provision for any other feedback loops, options 1 and 
5 are of limited value. 

Card 5 .- (((-feedback gains. FORMAT (6F10.4). This card may 
be left blank unless options 3 or 5 were selected. In this case, 
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six values of the gain K|,. must be specified, for wh\ch th,^ program 

4 > 

will determine the closed loop poles. This option, because of the 
limit to one feedback loop, is of little value. 

The user may analyze several flight conditions in a single 
computer run by stacking run data decks for several cases. Note 
that the run is terminated by the code word END in the title card. 


Output of Program ATDYN 

The output olf program ATDYN is best explained by using an ejr* 
ample of «n actual priritout of ATDYN, as ft is provided in Figure 28 

After printing out the initial conditions as specified by the 
input data to ATDYN, the aerodynamic coefficients as interpolated 
for this flight condition are printed. 

Next, the array of the coefficients are printed. (See sec- 
tions on linearized equation of motion for the definition of these 
coeffi cients . ) 


The options selected for running ATDYN are printed next (see 
description of input card 2). 

The "longitudinal matrix" is the matrix as specified by 

the first order matrix differential equation 

* = ''long * ®long 

with the X vector being ('u, 9, 0, a). This form assumes that the 
equations are known in first order form and Appendix A contains the 
formula for the coefficients of and as they are printed 

out here. 


Next, the three zeros (roots of the numerator polynomial) of 
the transfer function ol{s)/6^{s) are printed. These zeros were 
obtained by use of the subroutine ALFZER and POLFAC which requires 
an explicit first-order expression for the transfer function. The 
lengthy algebraic calculations to obtain this traasler function in 
its desired form are not included in this report^ 
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Figure 28 (Concluded) 
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The longitudinal roots are printed next in the form: 

Root Number Real Part Imaginary Part 

are the eigenvalues of the longitudinal matrix printed out above. 
These roots were obtained by the subroutine EISPACK and the error 
message following these roots, when zero, indicates that there were 
no numerical difficulties encountered by the EISPACK routines in 
finding the eigenvalues of the matrix. 

The reittli ndiar of the longitudinal dynamics printout was gener- 
ated by the subroutine TRANS. A detailed description of how to 
interpret the printout of TRANS is provided in the section des- 
cribing the program LATPOL. The matrices for subroutine TRANS were 
set up according to the linearized equations as listed in the sec- 
tion “Linearized Equations for the Longitudinal Direction" as fol- 
lows: 


■u 

— 
* a 

6 

«s 

sK^+Kg 

SK3+K4 

sKg+Kg 

^7 

00 

sKg+Kio 

sK,, + K,2 


hi 

sKi5+Ki6 

s^K|7+sKi8 

Kl g 


The numerator polynomial where column 2 is replaced by col- 
umn 4 is therefore the numerator polynomial of the transfer func- 
tion 'a(s)/6g(s). 

The denominator polynomial (for the pbles of the transfer 
functions) is obtained from the determinant of the matrix consis- 
ting of the three rows and the first three columns. (The printout 
annotates this with "Column 4 replaced by column 4"). 

Note that both the zeros and poles of the a(s)/6g(s) transfer 


134 



I 


function as calculated iby the first order differential equations 
and by the subroutine TRANS are identical. 


The program ATDYN then sets up the equations of motion for the 
1 ateral -di recti onal axes in first order differential equation form 

^ ^ ^lat ^ * ^lat,l ^ ^lat,2 


The x-vector is (4), 4> , 4), 4;, g). The derivation of the elements in 
the above matrices is too lengthy td be included in this report. 

The matrix A-j^^ and the vectors -j and 2 printed 

next, followed by the roots of the characteristic equation. 

Next, tile resoH? of usi subroutine TRANS are printed. The 
matrices were set up as follows (compare section “Linearized Equa- 
tions for Lateral Direction"): 




3 

«a 

^R 

s KjQ+sKg^ 

s^K32+sK33 

*^34 

*^35 

*^36 

s ^37+5X23 

s^K39+sK4o 

•^41 

K42 

K43 

SK44+K45 

s>^46 + ^47 

sI<48''’K49 

*^50 

^51 


The numerator zeros (column 1 replaced by column 4) are there- 
fore the zeros of the transfer function <j>(s)/6.(s) while the roots 

a 

of the characteristic equation are obtained by replacing column 4 
with column 4. 


Note again the perfect agreement in the values of the roots 
of the characteristic equations as obtained by EISPACK (as eigen- 
values of the lateral matrix) and by TRANS (as zeros of the poly- 
nomial of the characteristic equation). 

Finally, the program ATDYN prints out the results of the sub- 
routine LATPOL. Remember, that LATPOL was written before the sub- 
routine TRANS was available. First printed are the polynomial co- 
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efficients (polynomial of the characteristic equation). There is a 
slight difference in the coefficient of s® between TRANS and LATPOL 
which may be due to round-off errors. Also printed are the zeros 
of the transfer function 4>(s)/6g(s) (they agree exactly with those 
found by TRANS), and the zeros of the transfer function <|)(s)/6p(s) . 

In addition to the printed output of ATDYN, a deck of cards is 
punched containing the array of the coefficients K^. , i=l to 51. 

The format for this punched output is 8F10.6, so the entire deck 
for the K-array consists of nine cards. 


Program lONLQP 

(Closed Loop Longitudinal Axis) 

Purpose . - To obtain closed loop roots and zeros of the linear- 
ized longitudinal axis. 

Input .- K-array , Deck (K^ through ) 

7 cards FORMAT (8F10.6) 


Title Card: 

1 card 
Longitudinal Gain Card: 

This card contains one set of gains, that is 


FORMAT (78A1) 


'^P..* *"'1 » '^R ’ ^Lon 


a 


a 


a 


1 Card 

Lateral Gain Card: 


FORMAT (4F10.2) 


This card is read by LONLOP, but any information on it is ignored. 
(Should be blank card.) 


Run is terminated when title card contains END in columns 1, 

2 and 3. 

Output . - Printout of longitudinal closed loop roots and zeros 
as printed by subroutine TRANS and explained in writeup of ATDYN. 
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Method s- the program sets up the matrix of the closed lobp 
longitudinal control loop as explained below; 

For the purpose of stability analysis of the longitudinal con- 
trol loop these equations were combined with the control equations 
portrayed in FigurelS in a 6 x 7 matrix format suitable for analy- 
sis using the subroutine TRANS. This matrix takes the form 



1 

2 

3 

4 

5 

6 

7 


- 'u 

a 

0 

«s 

^sc 

«e 

“c 

1 

( K*| S+K2^ 



(-K 7 ) 

lo) 


(0) 

2 

(Kg) 

(KgS+Kio) 

(K^^s + K^2^ 

(-K13) 

(0) 

(0) 

(0) 

3 

(K,4) 


(S^^s^+K^gs) 

(-Kig) 

(0) 

(0) 

(0) 

i 

(0) 

(A) 

(0) 


c-n 

: CO) 

m 

■ 

(0) 

(0) 

(0) 

(0) 

(s) 


(0) 

6 

(0) 

(1 + Kp Kp s) 
a a 

(0) 

(0) 

(0) 

(1) 

(1) 


The variable is the commanded perturbation in angle of attack. 


The first run should usually have blank cards for data. The 
characteristic equation with zero input gains should contain the 
aircraft longitudinal poles just as they were printed in ATDYN, and 
this serves as a check run. 


Program LABEZR 

(Open Loop Zeros and Poles B(s)/5j^(s)) 

Purpose . - To calculate the open loop zeros and poles of 
0(s)/6j^(s) for one specific flight condition. 
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Input . * K-array deck: (K^ through Kg^ ) 

7 cards FORMAT (8F10.6) 

Title card: 

1 card FORMAT (78A1) 

Output . - Open loop zeros and poles of 8(s)/6p(s). 

Method . - The matrix as given in the description of subroutine 
TRANS is set up and then TRANS is called. 

Program LATDYN (Closed 3-Loop) 

Purpose . - To calculate poles and zeros of the lateral dynam- 
ics with only the 3-loop closed. This is used to vary Kp and Kp 

3 B 

until satisfactory dutch roll characteristics are obtained. 

Input .- K-array deck (K^ through 

7 cards FORMAT (8F10.6) 

Title card: 

1 card FORMAT (8F10.6) 

Longitudinal gains 

1 card blank card 

Lateral gain Kp and Kp 

FORMAT (30X, 2F10.2) 

The runs are terminated if the title card contains END in columns 
1 through 3. 

Output . - Closed loop poles and zeros of the transfer function 
<l>(s)/6g(s). 

Method . - The following matrix is set up for use by TRANS: 
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<p 


B ' :■ 

«a 

*^30^ '*’'^31^ 

*^32^ ''''^33^ 

'^34'^'^36'^P ■^'^36'^p/r ^ 

P P P 

'^35 

*^37^ ■'■’^38^ 

K39S^+K4qS 

K 43 +K 43 KP +K 43 KP Kp s 

P P P 

K 42 


K46S+K47 

K49^K5,K .(KjiKp K 

P P P 

*^50 


The first run should use K„ = Kp =0 for comparison of the 

e B 

poles with the poles generated by ATDYN. 


Program LATLOP (Closed Loop Lateral) 

Purpose . - To calculate the closed loop poles and zeros of the 
transfer function 4) ( s ) /cj)^ ( s ) . 

Input .- K-array deck. 

7 cards FORMAT (8F10.6) 

Title card 

1 card FORMAT (78A1) 

Blank card 

1 card blank 

Lateral gains Kp , Kj , Kp , Kp , Kp , T|_^^ 

(j) (j) (|) 3 B 

1 card FORMAT (6F10.2) 

The runs are terminated by END in columns 1, 2 and 3 of the 
title card. 

Output.- Poles and zeros of the transfer function 4 )(s )/4)^(s ) . 

Method . - The following matrix, reflecti ng the bl ock diagram 
represented in Figure 29 is set up, for analy,sis by the subroutine 
TRANS. 
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Program LALINI 

(Lateral Linear Integration) 

Purpose . - To generate time histories of the linear, lateral, 
closed loop system. The aircraft dynamics are linear, but the 
nonlinear position and rate limit of the ai 1 erons-spoi 1 ers are used. 

Input . - The K-array deck. 


7 cards 
Title card 

FORMAT 

(8F10.6) 

1 card 
Blank card 

FORMAT 

(78A1 ) 

1 card 

blank 


Lateral gains Kp , K. , 
% (|> 

1 card 

Initial Conditions 

Kr . Kp , 

% % 
FORMAT 

•^Rg’ "^Lat 
(6F10.2) 

ijj, ((), ij<, 4), 6, 

(all in deg. and deg/sec) 

1 card 

Run Parameters 

DT, Tg^^, NPRNT 

FORMAT 

(6F10.3) 


where 

DT = integration stepsize 
Tend " time 

NPRNT= integration steps between printout 

Output . - The following variables are printed: 

Time ^ 4> ij; ({» 6 6 6 6 6 6 

®p ^c ®RL ® 

(in degrees and degrees/second) 
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INPUT-OUTPUT UNIT ASSIGNMENTS 


The program AML75 requires the assignment of the following 
logical FORTRAN I/O units to physical units. 


FORTRAN I/O 



LOGICAL UNIT 

PHYSICAL UNIT 


3 

Card Reader 


7 

Card Reader 


4 

Line Printer 


8 

Line Printer 


6 

Scratch Tape 1 


9 

Scratch Tape 2 


Fortran units 6 and 9 are 

used to create a data set 

used for 

subsequent plotting of command 

sequences and responses. 

These 

units may be bypassed at LRC. 

For running AML75 in mode 

2 (AML 

versus canned trajectory), the 

additional I/O requirement 

exi sts : 


11 Magnetic Tape with recorded trajectory 

The program ATDYN requires the following assignments: 

3 Card Reader 

4 Line Printer 

5 Card Punch 

The individual small programs for the linear analysis require 
the following assignments: 

3 Card Reader 

4 Line Printer 
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CONCLUSIONS 


The work performed under this contract clearly indicates that 
it is possible to replace the performance model presently used in 
the AML program by an aircraft model of the same complexity as the 
one presently representing the aircraft flown by the human pilot 
on the DMS. 

Limited experience with runs of the new equations of motion 
replacing the old performance model indicated, surprisingly, that 
the aircraft driven by the control system achieved higher cell 
values, that is, it was tactically superior to the old version. 

The computational effort associated with the portion of the 
program performing the control system simulation is quite small 
compared to the computational effort to solve the nonlinear equa- 
tions of motion with all its associated interpolation of aerody- 
namic coefficients. 

While the longitudinal control system is of a relatively con- 
ventional design, some novel and interesting concepts have been 
incorporated into the lateral axis control system. Here, the com- 
bination of a bang-bang control system for large errors and/or 
error rates with a linear control system for small error and error 
rates represents a concept which deserves consideration for imple- 
mentation in a real fighter aircraft. 
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RECOMMENDATIONS FOR FURTHER STUDIES 


Implementation of the New 
Control System on the DMS 

The following discussion points out some of the potential 
problems which might arise when the new AML program is implemented 
to run in real time on the DMS and suggests possible approaches for 
their solution. 

Required memory storage may exceed available core storage if 
the AML program in its present form is combined with the DMS pro- 
gram. Required memory could be reduced drastically by the follow- 
ing modifications to the AML program: 1) Use tables defining the 

aircraft and atmosphere already existing in the DMS program. This 
would not only result in a reduction of core requirements, but it 
would also guarantee identical data for the OMS-driven and the AML- 
driven aircraft. 2) Use the same interpolation routines as the 
DMS pi*d^r|i|te||is!!obtain aerodyne 

A s««on#and probably more serious problem might arisie in the 
required of time to perform the calculations for the equa- 

tions of mofi oh for the target. Computational requirements for 
just one single pass through the full equations of motion are many 
times larger than they were for the old equations of motion. The 
calculation of all stability and control derivatives was not pres- 
ent in the Old version. Furthermore, the calculation of the three 
moments L, IIICAnd N was not required by the old equations of motion 

and the three linear forces F, F and F were considerably simpler 
, X y z 

in the old version due to the assumption of no sideslip and the 
neglect of transient terms. To save computer time, the program is 
already set up to calculate certain stability and control deriva- 
tives only every fourth integration step. It is important to 
recognize that in order to run on the DMS in real time, those 
cycles which require the most computation time are the limiting 
factors. In the past, the limiting cycles were those during which 
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the AML program executed a tactical decision process. If now during 
those cycles already requiring the total available amount of com- 
putation time, the additional demand for the new equations of motion 
is made, not enough computer time will be available during those 
cycles. Several possible remedies for alleviating this difficulty 
can be considered* First, the programming of the AML tactical de- 
cision process could be improved by carefully sequencing the calcu- 
lations tnsucft .a way that each variable is calculated only once 
during a cycle; Tri gonometric functions might" be calcul ated "by ' 
table look up (many expressions involving trigonometric functions 
do not require the ten or better digit accuracy obtained by the 
standard library routines). Subroutine calls with very long argu- 
ment lists could be replaced by more efficient calls, and several 
other programming "tricks" could be used, may be even replacing 
the FORTRAN coding of certain critical portions by more efficient 
routines written in assembly language. 

Similar procedures could be used to speed up the solution of 
the equatiqn#;pf motion. The first step here should be ti| ^^Stagger" 
the calculat4cm^^^ stability and control derivatives. Macb iuimber 
and altitud^ within , say, one quarter of a second cbah|e not Enough 
to warraht t%cal cul ati on of all the terms which are functions of 
these two independent variables. Updating stability and control 
derivatives, a few of them during each integration cycle and may be 
none of them during the critical cycles where a tactical decision 
is made wbuTd spread out more evenly to requirements for computa- 
tion time. , . , 

Should all the above proposed improvements still not make the 
AML program to execute fast enough, one might consider to use sim- 
plified equations of motion during the three or four time intervals 
during which the tactical decision is made. For instance, one 
might consider the aerodynamic moments as remaining constant over 
these few integration steps. Some analytical study could indicate 
how valid such an assumption would be. 
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All the above proposed measures for speeding up the execution 
time of the AML program did not assume any changes to the DMS pro- 
gram itself. But there is a good chance that additional time could 
be gained by improving the computational efficiency of the DMS 
program. 

The final and ultimate step would be to use a second computer, 
working in paraTTel; with the 6600 and to make this computer dedicated 
to solving the -equations of motion for the AML program. A mini- 
computer of the size of a PDP-11 might be powerful enough to per- 
form that task. Solving the equations of motion and simulating 
the control system on an analog computer interfaced to the CDC 6600 
certainly would satisfy the computational requirements and relative- 
ly few signals would have to be transferred back and forth from the 
6600. However, it is this author's personal opinion that such a 
hybrid solution should only be attempted if everything else fails. 

Improved Performance Specifications 

. ' ■ / it- .. 

The authors of this report are convinced that the most signif- 
icant proff^ess in air-to-air combat in the coming years will not 
be made by developing new and improved tactics (we are talking still 
of only one-on-one engagements) but in developing better airframe 
control system combinations, especially for high angle of attack 
flight. A recent RFQ issued by the Air Force Systems Command (Ref. 
10) entitled "Identification of Key Maneuver-Limiting Factors" 
clearly shows the trend to better understand the flight at high 
angles of attack. Associated with the factors degrading high angle 
of attack maneuverability is the problem of how to specify aircraft 
performance characteristics at high angle of attack. A key ques- 
tion here is. How valid is it to replace the aircraft by an equiva- 
lent second-order system under these conditions? The inherent 
strong nonlinearities affecting these flight conditions would ap- 
pear to invalidate such an approach. How are handling qualities 
of an aircraft under these high angle of attack flight conditions 
best described? If it should turn out that fewer than six degrees 
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of freedom will not adequately describe high angle of attack flight, 
Mil Spec 8785 might simply no longer be applicable. Can a flight 
control system really alleviate these problems to such a degree 
that the aircraft then again could be described by a second order 
system? 

It is quite obvious that a lot more questions exist than there 
are answers, and a new approach to study this important subject ap- 
pears to be needed. 

Improved Control Strategies 

The simultaneous change of angle of attack and bank angle re- 
ceived some limited consideration in this study and a relatively 
simple algorithm for angle of attack reduction during large roll 
maneuvers was developed, as described in the section of the longi- 
tudinal control for the nonlinear system. We have no proof that 
this algorithm is optimal, the only statement we can make is that 
it works. to work quite well, but considerably 

more effort titan jiist engineering judgment should be put forth to 
refine this 11 It is certainly a key factor in executing 

simultaneous ro;11 and pitch maneuvers in minimum time; 

Since all the necessary computer programs to study this ques- 
tion are now available, it is strongly recommended that in the de- 
sign of a similar control system for another aircraft than the F-4, 
this problem be studied thoroughly. 

Improvements to the AML Program 

The AML program was developed specifically for an aircraft 
represented by a performance model. Certain restrictions on the 
aircraft motion were imposed, one of the most important ones that 
the aircraft fly in a maneuverpl ane . While the concept of the 
maneuverpl ane is without doubt of great value for setting up and 
for evaluating tri al , maneuvers ; however, it seems that little 
benefit is obtained by postulating that the aircraft, between 
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tactical decisions, fly in a maneuverpl ane . 

Other problems arising in the AML program show up drastically 
when flying with the new control system. The most pronounced of 
these difficulties occurs during pullups or dives through 90 degrees. 
Some ambiguities exist if the velocity vector has not yet made the 
transition through 90 degrees but the aircraft body axis, due to the 
angle of attack, already passed the singularity at 90 degrees. 

It is quite likely that increased usage of the AML program, 
with the new control system may reveal additional areas where the 
AML program could be improved. 

Decision Science, Inc. 

San Diego, California 
5 February 1976 
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TABLE I 

LIST OF SUBROUTINES USED BY AML75 


1 . 

AERCF 

25. 

LACNTR 

2. 

AERF4 

26. 

LOCNTR 

3. 

AERFUN 

27. 

MEINTL 

4. 

AMLVS3 

28. 

NLPRNT 

5. 

GLOSS 

29. 

NONLIN 

6 . 

CMDSEQ 

30. 

NORPLN 

7. 

CMTRX 

31 . 

PAIRCR 

8. 

CSRHO 

32. 

PCELL 

9. 

CSRHOX 

33. 

PRCELA 

10. 

DEBGCT 

34. 

PRCELT 

11 . 

DIRCOS 

35. 

PRESR 

12. 

EQMOTA 

36. 

PRTF4 

13. 

EQMOT 

37. 

PRTNEW 

14. 

ERMSG 

38. 

REACTA 

15. 

EXTRA 

39. 

REACTT 

16. 

EXTRT 

40. 

RELGN 

17. 

FUHCl 

41 . 

RTSUB 

18. 

FUNC2A 

42. 

STATEA 

19. 

FUNC2D 

43. 

STATET 

20. 

FUNC2H 

44. 

THRTLA 

21 . 

FUNC2X 

45. 

TROTLT 

22. 

FUNC3 

46. 

TRYNXA 

23. 

GETCOM 

47. 

TRYNXT 

24. 

GETRXN 
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TABLE II 

LIST OF PROGRAMS AMD SUBROUTINES FOR 
LINEAR ANALYSIS 


1. Program ATDYN uses the following subroutines not yet listed 


2. 

ALFZER 

under AML75 

3. 

BALANC 


4. 

ELMHES 


5. 

FACTOR 


6. 

HQR 


7. 

KS 


8. 

LATPOL 


9. 

POLDRW 


10. 

POLFAC 


11. 

POLRT 


12. 

RG 


13. 

TRANS 



The following five main programs require no subroutines not listed 
earl i er : 

14. LONLOP 

15. LABEZR 

16. LATDYN 

17. LATLOP 

18. LALINI 
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APPENDIX A 


Transformations of the Equations of Motion 

Linear differential equations may be transformed in many dif- 
ferent ways for the purposes of control analysis. In the cases 
that follow the transformations made computation more efficient 
and ens^fed finding the numeirator dynamics or zeros of the transfer 
function of interest. 

dif ferei|t1 al equiti ons of motion for the longitudinal or 
pitch equailOns (l)r contain first and second derivatives of 

the fo,|&*r' ytriablo®^ For numerical integration and state 

variable methods of analysis it is desirable to uncouple the deriv- 
atives resulti ng in one differential equation of first order for 
each variable present. 

The first-order differential equation form of the equations 
pitch of motion may be obtained by algebraic elimination of coupled 
derivative terms from the linearized pitch dynamics, equations ( 1 ) . 
If a state vector 

'u 

^2 - 6 

X 3 = 0 

X4 = a 

is defined, equations (1) assume the form 

KiX-^ ^ ^2^1 ^3^4 ^ ^ ^ ” ^7^s (x— axis) 

KgX-j + KgX^ + ^1 Q ^4 ^ 11^3 ^12^2 ~ ^13^s (z-axis) 


K14X1 + 


^ 15^4 ^ 


K16X4 + 


K17X3 + K|gX3 = 


(pi tch) 


A1 



The z-axis equation may be solved for and the result used 
to eliminate x^ from the x-axis and pitch equations. The resulting 
equations are: 




+ K V V (K7-K-J ^ )6 


1 


4 ’'3 Kg '‘"7 ‘'3 Kg ' “s 


Xg = X 3 


1 


K 


17 


^15 ^12 
'17 “^9 


K, 


^3 ^ K,-,^‘‘^14'^'^15 * K,'! k1^ ^^2 ^ S'^'^l 5 1^^^3 


'17 


K, 


• ^8 ^12 ^11 ^10 ^13 

"V 'V ‘V *4 ^ 1:;“ ®s 


and 


These equations define the coefficients of the matrix 
Biong in the equation 


* ' Nong + “^long «s 

used in ATDYN and LONLOP. 
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